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Abstract 

Solid-state  quantum  sensors  are  attracting  wide  interest  because  of  their  sensitivity 
at  room  temperature.  The  spin  properties  of  individual  nitrogen  vacancy  (NV)  color 
centers  in  diamond  make  them  outstanding  nanoscale  sensors  of  magnetic  fields,  elec¬ 
tric  fields,  and  temperature  under  ambient  conditions.  Using  large  ensembles  of  NV 
centers  leads  to  increases  in  sensitivity  with  reduced  spatial  resolution.  A  power- 
efficient  device  for  the  implementation  of  diamond-based  sensors  using  large  ensem¬ 
bles  (>  1010)  of  NV  centers  with  sub-millimeter-scale  spatial  resolution  is  described 
and  demonstrated.  The  challenges  of  working  with  large  ensembles  of  defects  are 
be  discussed,  as  well  as  several  strategies  for  combatting  sensitivity-limiting  ambi¬ 
ent  temperature  fluctuations.  Magnetometry,  temperature  sensing,  and  timekeeping 
applications  are  explored. 
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rection  of  applied  magnetic  held,  where  B\\  =  Bz  and  B±  =  sjB\  +  B2y. 
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1-2  ODMR  (a)  shows  a  single  NV  resonance  with  no  applied  magnetic 
held  and  at  room  temperature  (To)-  The  resonance  linewidth  (5f)  and 
normalized  fluorescence  difference  between  the  bright  and  dark  states 
( C )  are  indicated.  The  addition  of  a  magnetic  held  aligned  to  the  NV 
axis  causes  the  resonance  to  split  (b)  while  a  change  in  temperature 
causes  a  global  shift  (c).  This  behavior  can  be  generalized  to  the  four 
sub-ensembles:  the  splitting  of  each  set  of  resonances  is  proportional  to 
the  projection  of  the  applied  magnetic  held  onto  that  sub-ensemble’s 
axis .  24 
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2-1  Light  Trapping  Diamond  Waveguide  (a)  Standard  single-pass 
laser  excitation  exhibiting  a  short  path  length  of  interaction  with  NV 
centers  as  compared  to  the  LTDW.  (b)  For  an  equivalent  excitation 
power,  the  LTDW  increases  the  path  length  of  the  excitation  beam  by 
up  to  several  orders  of  magnitude,  (c)  The  colour  CCD  image  of  the 
LTDW,  excited  by  a  green  pump  laser  of  70  mW,  shows  bright  fluores¬ 
cence  without  any  spectral  filtering,  (d)  Simulation  of  the  beam  path 
in  the  LTDW  with  a  =  500  pm  input  facet,  3  mm  side  length,  100  pm 
diameter  Gaussian  profile  input  incident  at  11.6°,  and  an  absorption 
constant  of  a  =  0.45  cm-1 .  33 

2-2  Optical  properties  of  the  LTDW  (a)  Pump  photon  flux  is  plot¬ 
ted  against  detected  photon  flux  for  the  LTDW  (blue)  and  for  a  single 
pass  through  the  long  side  of  the  device  (red).  Error  bars  are  con¬ 
tained  within  the  data  points,  (b)  Maximum  optical  path  length  (un¬ 
constrained  by  absorption)  is  plotted  as  a  function  of  input  angle  into 
the  LTDW  structure  with  a  150  pm  input  facet  and  standard  beam  di¬ 
vergence  in  blue.  The  excitation  laser  intersects  the  input  facet  50  pm 
from  the  center,  resulting  in  asymmetry  in  the  plot.  For  comparison, 
the  maximum  optical  path  length  achievable  in  a  single-pass  config¬ 
uration  through  the  top  of  the  sample  is  also  plotted  in  black,  and 
through  the  side  of  the  sample  is  plotted  in  red.  In  samples  with  suffi¬ 
ciently  high  NV  density,  the  excitation  beam  would  be  absorbed  before 
achieving  the  maximum  path  length  allowed  by  the  LTDW  structure.  34 

2-3  (a)  A  MATLAB  simulation  example  of  a  pattern  where,  after  32  re¬ 

flections,  the  beam  exits  the  sample  at  an  angle  adequately  separated 
from  the  input  path  to  enable  transmission  measurements,  (b)  An 
example  of  a  path  that  spatially  covers  the  majority  of  the  diamond 
sample,  (c)  A  three-dimensional  rendering  of  the  LTDW  from  the  Ze- 
max  simulation  featuring  a  diverging  Gaussian  beam.  Beam  divergence 
is  artificially  exaggerated  by  two  orders  of  magnitude  for  clarity.  ...  35 
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2-4  Atomic  force  microscopy  characterization  of  surface  polishing 


42 


3-1  ODMR  spectra  and  level  structure  (a)  CW  electron  spin  reso¬ 
nance  with  external  magnetic  held  aligned  along  a  single  (111)  crystal 
axis.  The  inset  depicts  the  diamond  lattice,  where  the  blue  circle  rep¬ 
resents  the  nitrogen  atom  and  the  green  circles  show  the  four  possible 
orientations  for  the  vacancy  in  the  tetrahedral  crystal  lattice,  constitut¬ 
ing  the  four  sub-ensembles  of  NVs.  (b)  Lock-in  output  corresponding 
to  signal  in  (a).  The  scale  factor  (V/Hz)  is  provided  by  a  linear  fit 
around  each  of  the  two  red  points.  Temperature  and  magnetic  held 
shifts  are  measured  independently  by  tracking  both  rns  =  ±1.  (c) 

Energy-level  diagram  of  diamond  NV  center  showing  radiative  (solid 
lines)  and  non-radiative  (dotted  lines)  transitions .  46 


3-2  Sensitivity  measurements  (a)  Noise  density  of  the  in-phase  lock-in 
ampliher  signal  at  the  cc+  resonance  (see  inset)  converted  to  units  of 
magnetic  held  for  a  time  constant  of  100  ms;  low  frequency  noise  is 
attributed  to  thermal  effects.  (b)ESeparation  of  magnetic  held  and 
temperature  effects.  Green  and  blue  lines  show  the  sum  and  difference 
of  the  intensity  signals  from  lower  and  higher  frequency  resonances  of 
the  (111)  orientation  with  a  time  constant  of  10  ms,  decoupling  the 
temperature  and  magnetic  held  drifts  over  100  s,  respectively.  A  20-s- 
period  square  wave  of  amplitude  190  nT  is  applied  with  external  coils 
and  is  clearly  visible  in  the  blue  trace.  The  noise  on  a  single  10  s  step 
(0.93  nT)  extrapolates  to  ~  290  pT  Hz'1/2  when  accounting  for  the 
10  Hz  equivalent  noise  bandwidth  of  the  lock-in  ampliher.  The  green 
trace  shows  drift  attributed  to  environmental  effects;  aliasing  effects 
have  been  removed  (less  than  200  ms) .  47 
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3-3  Comparison  of  the  LTDW  (dark  blue  overlay)  with  other  mag- 
netometry  technologies  Particularly  with  regard  to  target  applica¬ 
tions  such  as  geomagnetics  [1]  shown  in  green  (corresponding  to  the 
right  axis),  the  LTDW  achieves  improved  sensitivity  with  respect  to 
other  diamond-based  sensors,  shown  in  bine  (corresponding  to  the  left 
axis),  in  the  important  low- frequency  band  [2],  The  LTDW’s  perfor¬ 
mance  is  competitive  with  other  magnetometer  platforms  [3] ,  shown  in 
light  bine .  52 


4-1  (a)  Axial  and  transverse  magnetic  held  directions  with  respect  to  NV 

axis.  As  B  j_  is  increased  from  35  G  (b)  to  45  G  (c),  the  rrij  =  0  energy 
levels  are  separated  from  the  rrij  =  ±1  energy  levels,  removing  the 
degeneracy  which  results  in  an  electron-nnclear  spin  hip  anti-crossing, 
which  is  detailed  in  the  inset  of  (b).  Non-dressed  states  are  indicated 
with  colored  lines  and  dressed  states  are  indicated  with  black  lines.  .  57 

4-2  (a)  Anti-crossing  coupling  strength  with  increasing  B±.  (b)  Magnitude 

of  small  axial  B  helds  at  which  the  centers  of  these  anti-crossings  are 
observed.  In  both  plots,  the  double-electron  spin  hip  anti-crossing  is 
marked  E  and  the  electron-nuclear  spin  hip  anti-crossing  is  marked  N. 

At  helds  greater  than  B±  ~  40  G,  anti-crossing  E  continues  to  increase 
quadratically  and  anti-crossing  N  goes  to  zero  as  the  energy  levels  no 
longer  overlap .  60 

4-3  (a)  Experimental  setup,  (b)  The  four  orientation  sub-ensembles  of 

NV  centers  are  shown  with  an  applied  magnetic  held,  (c)  An  anti¬ 
crossing  strength  of  ~  2.6  MHz  is  seen  in  sub-ensemble  1  when  the 
perpendicular  component  of  the  B  held  (B i;j_)  is  equal  to  30  G  and  the 
axial  component  of  the  B  held  (Him)  is  varied  around  zero  held.  Inset 
shows  a  high  resolution  spectra  of  resonance  1  and  1' .  61 
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4-4  Plots  of  transition  energies  as  a  function  of  axial  magnetic  field  under 
fixed  transverse  magnetic  field.  Transitions  for  B±_  =  35  G  and  B±  = 

45  G  are  shown,  which  correspond  to  the  energy  levels  in  Fig.  4-1. 
Greyscale  indicates  normalized  fluorescence.  Double-electron  spin  flip 
anti- crossings  of  coupling  strength  Eg  are  seen  at  Ep  =  0,  ±A\\h/geHB  — 

±0.75  G.  ESR  in  diamond  with  a  natural  abundance  of  carbon  isotopes 
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Chapter  1 


Introduction  and  Background 


Solid-state  quantum  sensors  are  attracting  wide  interest  because  of  their  excellent  sen¬ 
sitivity  at  room  temperature.  In  particular,  the  spin  properties  of  individual  nitrogen 
vacancy  (NV)  colour  centers  in  diamond  [4,  5,  6]  make  them  outstanding  sensors  of 
magnetic  fields  [7,  8,  9,  10,  11,  12,  13,  14,  15],  electric  fields  [16],  strain  [17],  and 
temperature  [18,  19,  20]  under  ambient  conditions.  By  employing  an  ensemble  of  N 
unentangled  NVs,  a  factor  of  up  to  a /N  improvement  in  sensitivity  is  available  [14,  21]. 

This  thesis  work  began  in  pursuit  of  a  diamond  electron-spin  clock  for  space 
navigation  and  communication,  as  a  NASA  Space  Technology  Research  Fellowship. 
The  goal  was  to  implement  a  theoretical  proposal  for  timekeeping  in  diamond  [22], 
For  the  clock  stability  to  be  competitive  with  modern  compact  atomic  clocks  (~ 
10 -11r1/2)[23],  an  extremely  large  ensemble  of  NVs  is  required.  Furthermore,  dia¬ 
mond’s  sensitivity  to  temperature  fluctuations  presents  a  limitation  on  the  stability 
of  the  diamond-based  clock. 

New  techniques  were  required  to  efficiently  excite  a  large  NV  ensemble  and  to 
collect  the  optical  signal.  In  this  thesis,  after  a  brief  introduction  (chapter  1),  a 
light-trapping  diamond  waveguide  (LTDW)  geometry  is  introduced,  whose  excep¬ 
tional  excitation  efficiency  and  good  signal  collection  enables  in  excess  of  5%  conver¬ 
sion  efficiency  of  pump  photons  into  optically  detected  magnetic  resonance  (ODMR) 
[24]  fluorescence,  a  three  orders  of  magnitude  improvement  over  previous  single-pass 
geometries  (chapter  2).  This  dramatic  enhancement  of  ODMR  signal  enables  pre- 
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cision  broadband  separable  measurements  of  magnetic  field  and  temperature  in  the 
low-frequency  range  inaccessible  by  dynamical  decoupling  techniques  (chapters  3). 
It  enables  the  description  of  a  newfound  level  anti-crossing  in  the  diamond  energy 
level  structure  (chapter  4).  Preliminary  work,  using  this  device  for  pulsed  operation 
and  progress  towards  a  diamond-based  clock  is  detailed  in  chapter  5,  including  the 
continued  limitations  and  several  methods  for  overcoming  the  always-present  ambient 
temperature  fluctuations.  Finally,  an  outlook  for  the  held  of  ensemble-based  diamond 
sensing  and  timekeeping  is  discussed  in  chapter  6. 


1.1  Diamond  nitrogen  vacancy  defect  centers 

1.1.1  Description  of  the  physical  system 

Diamond  is  composed  of  a  tetrahedral  lattice  of  covalently-bonded  carbon  atoms  (Fig. 
1-la).  Defects  in  diamond  materials,  both  those  mined  from  the  earth  and  those 
created  in  labs,  have  been  studied  for  decades  [25]  and  dozens  of  known  defects  have 
been  characterized  [26].  Defect  centers  in  diamond  are  often  called  color-centers,  since 
they  can  give  this  normally  transparent,  colorless  material  a  visible  hue,  depending  on 
the  composition  of  the  defect.  Diamond  is  almost  98%  spin-neutral  12C  which  creates 
a  stable  vacuum-like  environment  surrounding  the  defect.  This,  through  engineers 
growth  with  isotopically  pure  12C,  can  be  improved  to  minimize  spin-bath  noise.  The 
NV  center  consists  of  a  carbon  atom  replaced  by  a  nitrogen  atom,  adjacent  to  a 
missing  carbon  atom  (Fig.  1-lb).  This  can  be  regarded  as  a  “trapped  ion”  in  the 
crystal  lattice.  There  are  two  optically  active  charge  states:  NV°  and  NV-,  but 
for  the  purpose  of  this  work  the  term  NV  refers  to  the  NV-  defect.  Additionally, 
in  diamond’s  tetrahedral  lattice,  there  are  four  possible  orientations  in  which  NV 
defects  can  be  positioned,  creating  four  sub-ensembles  of  defects.  These  four  stable 
sub-ensembles,  each  109.5°  separated  from  the  next,  enable  vector  measurements  of 
applied  fields. 

The  electronic  spin  state  of  the  NV  center  can  be  optically  polarized  using  a 
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Figure  1-1:  NV  structure  (a)  shows  the  defect-free  C-C  tetrahedral  covalent  bond 
structure,  (b)  shows  the  C-N  bond  with  the  vacancy.  The  lighter  colored  electron  is 
captured  from  the  lattice,  giving  the  NV  a  negative  charge,  (c)  shows  the  geometry  of 
the  NV  center,  indicating  the  the  direction  of  applied  magnetic  held,  where  B\\  =  Bz 
and  B±  =  sjB\  +  B^.  (d)  outlines  the  relevant  energy  levels  for  performing  ODMR. 


green  laser  and  read  out  by  monitoring  its  fluorescence  while  sweeping  an  applied  RF 
source.  Both  of  these  characteristics  are  attractive  properties  for  an  optical  sensor. 
When  driven  on  resonance,  the  fluorescence  level  decreases.  The  exact  frequencies 
of  these  resonances  is  determined  by  the  sensor’s  environment,  and  recording  these 
resonances  shifting  gives  insight  into  magnetic  held  and  temperature  experienced  by 
the  sensor. 


1.1.2  ODMR  in  diamond 

Further  insight  into  the  nature  of  the  observed  frequency  shifts  can  be  obtained  by 
modeling  the  energy  of  the  ground  state  triplet  sub-levels  in  the  weak  held,  low  strain 
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approximation  with  the  following  Hamiltonian: 


7-Lgg  —  DgsSl  +  gsfiBBzSz  + 


(1.1) 


where  Dgs  is  the  ground  state  crystal  held  splitting,  gs  is  the  Lande  g-factor,  gB  is  the 
Bohr  magneton,  Bz  is  the  component  of  the  magnetic  held  along  the  NV  axis  (Fig. 
1-lc) ,  A\\  is  the  hyperhne  tensor  describing  interactions  with  local  nuclear  spins,  and 
Sz  is  the  spin  projection  onto  the  z  axis  with  eigenvalues  0,  ±h,  corresponding  to 
ms  =  0,±1  [27].  Figure  1-ld  shows  the  relevant  sub-levels. 
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Figure  1-2:  ODMR  (a)  shows  a  single  NV  resonance  with  no  applied  magnetic  held 
and  at  room  temperature  (T0).  The  resonance  linewidth  ( Sf )  and  normalized  fluores¬ 
cence  difference  between  the  bright  and  dark  states  ( C )  are  indicated.  The  addition 
of  a  magnetic  held  aligned  to  the  NV  axis  causes  the  resonance  to  split  (b)  while  a 
change  in  temperature  causes  a  global  shift  (c).  This  behavior  can  be  generalized  to 
the  four  sub-ensembles:  the  splitting  of  each  set  of  resonances  is  proportional  to  the 
projection  of  the  applied  magnetic  held  onto  that  sub-ensemble’s  axis. 


The  triplet  ground  state  consists  of  a  degenerate  ms  =  ±1  state  separated  from 
the  zero  state  by  Dgs.  Quantities  such  as  temperature,  crystal  strain,  and  external 
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electric  fields,  cause  shifts  in  Dgs.  We  assume  that  the  intrinsic  strain  across  the 
sample  does  not  change  over  the  course  of  the  measurement  and  that  observed  shifts 
in  the  zero  field  splitting  (ZFS)  are  caused  by  temperature  changes. 

532  nm  green  light  is  used  to  excite  the  NV  from  the  ground  state  to  the  excited 
state,  which  fluoresces  in  the  red  as  it  returns  to  the  ground  state.  NVs  in  the  ms  =  ±1 
state  can  follow  a  non-spin  conserving,  non- radiative  path  to  a  shelving  state  where  an 
IR  photon  is  emitted  and  the  NV  returns  to  the  zero  state.  In  this  way,  green  light  can 
be  used  to  polarize  the  NV  into  the  zero  state.  Microwave  excitation  on  resonance  at 
Dgs  drives  NVs  from  the  ms  =  0  state  into  the  ms  =  ±1  state,  resulting  in  a  decrease  in 
observed  fluorescence  (Fig.  l-2a).  Axial  magnetic  fields,  which  appear  in  the  Sz  term, 
lift  the  ms  =  ±1  degeneracy,  allowing  the  effects  of  magnetic  fields  and  temperature 
to  be  separated  to  first  order  when  monitoring  both  halves  of  the  resonance.  The 
difference  between  the  two  resonances  is  proportional  to  the  applied  axial  magnetic 
field  (Fig.  l-2b)  and  the  change  in  average  of  the  two  resonances  is  proportional  to 
the  change  in  temperature  (Fig.  l-2c).  This  spin-dependent  fluorescence  is  the  basis 
for  using  the  diamond  NV  as  a  sensor.  When  the  device  is  continuously  driven,  it’s 
referred  to  as  continuous  wave  (CW)  ODMR.  The  laser  and  RF  excitation  can  also 
be  pulsed  for  improved  signal  to  noise  (C')and  decreased  linewidth  (<5/). 


1.2  Brief  overview  of  previous  work  in  diamond-based 
sensing 

1.2.1  Magnetometry 

Commonly  used  commercial  technologies  like  flnxgate  magnetometers  and  GMR- 
based  systems  require  careful  periodic  calibration,  dne  to  environmental  and  material 
changes  in  the  sensors.  This  maintenance  is  impractical  for  remote  operation.  Con¬ 
versely,  diamond-based  magnetometry  is  based  on  a  splitting  between  resonances  that 
is  entirely  dependent  on  fundamental  physical  quantities,  which  eliminates  the  need 
for  calibration.  Furthermore,  this  splitting  can  be  separated  from  the  global  reso- 
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nance  shift  caused  by  temperature  changes,  making  a  diamond-based  sensor  much 
more  reliable.  A  diamond-based  solid-state  magnetometer  also  has  the  potential  to 
be  a  much  more  compact,  robust  system  than  similar  atomic  magnetometers  as  no 
vacuum  system  is  required. 

Existing  diamond-based  magnetometry  systems  can  be  broken  into  two  main  cat¬ 
egories:  low  and  high  frequency  operation,  where  the  border  is  loosely  defined  to  be 
around  100  Hz.  To  date,  sensitivity  records  in  diamond-based  magnetometry  have 
been  set  in  the  high  frequency  regime  [28] ,  taking  advantage  of  dynamical  decoupling 
(pulsed)  techniques  to  get  away  from  low  frequency  environmental  noise.  However, 
the  majority  of  magnetometry  applications  exist  in  the  lower  frequency  regime,  in¬ 
cluding  geomagnetics  [1]  and  medical  applications  [29,  2,  3].  Unlike  high-frequency 
magnetometers,  which  are  exclusively  based  on  pulsed  operation,  low-frequency  mag¬ 
netometers  can  also  be  based  on  CW  ODMR,  as  described  in  section  1.1.2  where  the 
laser  and  RF  excitation  are  applied  continuously,  as  well  as  those  based  on  pulsed 
application  of  laser  and  RF  excitation.  Pulsed  operation  adds  complexity  but  has 
the  potential  to  result  in  increased  sensitivity  by  taking  advantage  of  the  NVs’  long 
coherence  time.  However,  the  use  of  lock-in  amplification  to  increase  signal  to  noise  in 
CW  ODMR  measurements  results  in  competitive  sensitivity  as  compared  with  pulsed 
operation. 

Recently,  defects  in  silicon  carbide  are  also  emerging  as  a  promising  material  for 
solid-state  magnetometry  [30,  31,  32],  Results  described  in  this  thesis  could  poten¬ 
tially  be  extended  to  work  in  bulk  SiC  devices.  However,  at  this  date,  long  spin 
coherence  times  at  room  temperature  of  these  defects  have  yet  to  be  demonstrated. 

1.2.2  Temperature  Sensing 

Unlike  diamond-based  magnetometry,  which  is  fully  based  on  fundamental  constants, 
diamond-based  temperature  sensing  relies  on  the  phenomenological  temperature  de¬ 
pendence  At  room  temperature,  this  is  found  to  be  ~  —75  kHz/K  [33].  Previous 
work  has  included  characterizing  a  polynomial  describing  how  this  quantity  changes 
with  respect  to  temperature  [34],  as  well  as  investigations  into  the  physical  mecha- 
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nism  behind  this  shift,  which  is  believed  to  be  a  combination  between  the  thermal 
expansion  of  the  diamond  and  electron-phonon  interactions  [35].  Pulsed  protocols 
designed  to  be  insensitive  to  magnetic  fields  [22]  have  also  been  used  to  measure 
temperature  [18,  19,  20]. 

While  there  are  applications  for  stand-alone  temperature  sensing,  the  main  mo¬ 
tivation  for  studying  the  temperature  dependence  in  NVs  is  to  remove  the  effects  of 
ambient  temperature  changes  from  room-temperature,  ambient-pressure  NV  magne¬ 
tometers. 


1.3  Working  with  ensembles  of  NVs 

1.3.1  Motivations 

The  sensitivity  of  a  measurement  of  a  CW  ODMR  line  shift  is  limited  by  a  factor 
V'esr  ^  cTTr  where  V  is  the  resonance  linewidth,  C  is  the  normalized  fluorescence 
difference  between  the  bright  and  dark  states,  and  R  is  the  photon  collection  rate. 
Linewidth  can  be  improved  to  some  extent  with  better  isolated  and  higher  quality 
NVs,  as  well  as  by  operating  at  low  temperatures.  This  can  be  one  to  two  orders 
of  magnitude  improvement,  at  the  expense  of  NV  density.  Normalized  fluorescence 
difference  can  be  improved  by  roughly  an  order  of  magnitude  by  optimizing  the  driv¬ 
ing  RF  and  laser  power  applied.  However,  by  increasing  the  number  of  NV  centers 
used  in  a  measurement,  the  R  value  can  be  increased  substantially,  leading  to  im¬ 
proved  sensitivity.  An  ensemble  of  1010  NVs  leads  to  a  improvement  of  five  orders 
of  magnitude.  This  sensitivity  advantage,  combined  with  some  of  the  advantages  of 
diamond-based  magnetometry  described  in  section  1.2.1  are  appealing  for  work  to¬ 
wards  a  mass-produced,  deployable  sensor  that  could  be  used  in  remote  and  isolated 
environments  without  the  need  for  periodic  recalibration. 

In  addition,  the  increased  sensitivity  achievable  through  ensemble-based  magnetic 
held  sensing  lends  itself  to  many  existing,  mature  applications,  where  single-NV  mag¬ 
netic  held  sensing  remains  limited  to  the  lab  setting.  In  particular,  geoscience  and  ob- 
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ject  detection  are  two  application  spaces  whereas  ensemble-based  sensing  shows  great 
potential.  Others  include  biomagnetism,  bioassays,  condensed  matter  physics,  mate¬ 
rial  science,  biothermometry,  bulk  magnetometry  for  surveying,  and  hyper-polarized 
media  for  NMR. 

1.3.2  Challenges 

The  challenges  facing  ensemble-based  NV  sensing  can  be  broken  into  two  categories: 
material  properties  and  experimental  design.  As  more  NVs  are  introduced  into  a 
diamond  sample,  either  during  the  growth  process  or  implanted  later,  there  is  an 
increase  lattice  damage  and  unwanted  defects.  This  leads  to  a  decrease  in  the  coher¬ 
ence  time  of  the  NV  which  negatively  affects  the  potential  sensitivity.  Unlike  a  single 
NV,  which  provides  atomic-scale  resolution,  the  resolution  provided  by  an  ensemble 
of  NVs  is  determined  by  the  volume  of  diamond  encompassing  the  ensemble.  This 
can  be  minimized  by  increasing  the  NV  density,  however,  this  comes  at  the  expense 
of  decreased  NV  quality.  The  devices  described  later  in  this  thesis  use  a  3  mm  x  3 
mm  x  0.3  mm  diamond  with  ~  0.1  ppm  NV  density.  Work  done  in  characterizing  this 
sample  will  be  discussed  in  chapter  2. 

A  basic  NV  ODMR  setup  requires  an  excitation  laser  and  RF  input  with  adequate 
power,  as  well  as  a  way  to  collect  the  signal-containing  fluorescence.  For  a  single  NV, 
a  common  configuration  is  to  have  the  excitation  laser  focused  onto  the  NV  and  a 
nearby  wire  or  antenna  providing  the  RF  excitation  -  uniformity  of  the  laser  and  RF 
excitation  is  not  a  concern.  However,  as  the  volume  of  diamond  that  the  NV  ensemble 
occupies  increases,  uniform  laser  and  RF  excitation  becomes  important,  as  both  affect 
the  contrast  of  the  signal,  particularly  when  working  below  saturation  [36] .  For  pulsed 
measurements,  this  becomes  even  more  important,  as  the  Rabi  frequencies  will  vary 
as  well.  The  light-trapping  diamond  waveguide  (LTDW)  geometry  that  is  described 
in  chapter  2  and  used  in  chapters  3,  4,  and  5  is  a  step  towards  more  uniform  laser 
excitation,  and  several  antenna  designs  are  described  that  work  towards  an  increase  in 
RF  excitation  uniformity.  Similarly,  imperfect  uniformity  in  an  applied  bias  magnetic 
held  leads  to  an  increase  in  resonance  line  width  and  a  loss  of  sensitivity.  Later 
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chapters  will  discuss  the  use  of  Halbach  arrays  and  Helmholtz  coils  to  address  bias 
held  uniformity. 

Another  challenge  in  working  with  NV-based  sensors  is  the  signal  collection  effi¬ 
ciency.  Diamond’s  high  index  of  refraction  traps  the  majority  of  the  fluorescence  via 
total  internal  reflection  (TIR).  Collection  through  the  same  path  as  excitation,  with 
the  addition  of  a  dichroic  filter  to  separate  the  excitation  and  signal  beams  based 
on  wavelength,  has  a  notoriously  low  efficiency,  often  less  than  1%.  Work  towards 
this  goal  with  single  and  small  ensembles  of  NVs  includes  efforts  using  solid  immer¬ 
sion  lenses  [37]  and  fiber-coupled  designs  [38].  In  ensemble-based  systems,  efforts 
include  early  work  with  diamond  waveguides,  placing  detectors  directly  adjacent  to 
the  diamond  [14]  and  the  use  of  parabolic  concentrator  lenses  [28]. 
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Chapter  2 


Light-trapping  Diamond  Waveguide 
Platform 


Portions  of  this  work  appeared  in  print  here:  [39] 

2.1  Introduction 

Solid-state  quantum  sensors  are  attracting  wide  interest  because  of  their  excellent 
sensitivity  at  room  temperature.  In  particular,  the  spin  properties  of  individual  ni¬ 
trogen  vacancy  (NV)  colour  centers  in  diamond  [4,  5,  6]  make  them  outstanding 
nanoscale  sensors  of  magnetic  fields  [7,  8,  9,  10,  11,  12],  electric  fields  [16],  and  tem¬ 
perature  [18,  19,  20]  under  ambient  conditions.  Recent  work  on  NV  ensemble-based 
magnetometers  [13,  14,  15],  inertial  sensors  [40],  and  clocks  [22]  has  employed  N 
unentangled  colour  centers  to  realize  a  factor  of  up  to  \/rN  improvement  in  sensitiv¬ 
ity  [14,  21].  However,  to  achieve  this  potential  sensitivity  enhancement,  new  tech¬ 
niques  are  required  to  excite  efficiently  and  to  collect  the  optical  signal  from  large 
NV  ensembles.  Here,  we  introduce  a  light-trapping  diamond  waveguide  (LTDW)  ge¬ 
ometry  whose  exceptional  excitation  efficiency  and  good  signal  collection  enables  in 
excess  of  5%  conversion  efficiency  of  pump  photons  into  optically  detected  magnetic 
resonance  (ODMR)  [24]  fluorescence,  a  three  orders  of  magnitude  improvement  over 
previous  single-pass  geometries.  This  dramatic  enhancement  of  ODMR  signal  en- 


31 


ables  precision  broadband  measurements  of  magnetic  field  and  temperature  in  the 
low-frequency  range  inaccessible  by  dynamical  decoupling  techniques. 

The  NV’s  low  absorption  cross  section  [41]  has  thus  far  prevented  high  conver¬ 
sion  efficiency  from  excitation  power  to  NV  ODMR  signal  in  bulk  diamond  samples 
using  the  typical  single-pass  excitation  geometry  illustrated  in  Fig.  2-la.  For  NV 
densities  around  1015  cm-3  [8,  42],  for  which  the  NV  spacing  in  the  diamond  lattice 
is  sufficiently  sparse  to  maintain  long  spin  coherence  times  [9],  up  to  meter- long  path 
length  would  be  necessary  for  significant  excitation  absorption  (see  section  2.2.2).  Ef¬ 
ficient  collection  has  been  demonstrated  from  NVs  in  single-pass  excitation  geometries 
[13,  14],  but  absorption  is  less  than  1%  for  a  typical  300  pm  thick,  electron-irradiated, 
type  Ha  diamond  sample.  Wide-held  microscopy  implemented  with  charge-coupled 
device  (CCD)  cameras  can  collect  fluorescence  from  ensembles  of  up  to  103  NV  cen¬ 
ters  [43 1  in  the  focal  volume  of  the  objective,  but  suffers  from  low  collection  efficiency 
with  the  majority  of  fluorescence  emission  trapped  within  the  diamond  due  to  total 
internal  reflection  (T1R)  at  the  diamond-air  interface.  External  Fabry-Perot  cavities 
could  increase  the  optical  depth  [15]  for  green  excitation,  but  would  introduce  band¬ 
width  restrictions  as  well  as  the  additional  complexity  of  stabilized,  narrow- linewidth 
excitation  lasers. 


2.2  Device  Design 

2.2.1  Description  of  the  LTDW 

We  have  designed  and  realized  a  new  ensemble  measurement  scheme  that  overcomes 
these  limitations.  The  LTDW  consists  of  a  rectangular  diamond  slab  with  a  small 
angled  facet  at  one  corner  for  input-coupling  of  the  pump  beam,  as  illustrated  in 
Fig.  2-lb.  The  input  facet  has  a  length  a  at  an  angle  45°  relative  to  the  square 
sample  sides,  allowing  pump  light  to  couple  into  the  structure  while  being  confined 
by  T1R  (6  >  6C  =  24.6°)  on  the  other  surfaces.  Fig.  2-lc  shows  an  optical  image  of 
a  pump  beam  coupled  at  an  11.6°  angle  into  the  LTDW  (3  mm  x  3  mm  x  300  pm). 
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Figure  2-1:  Light  Trapping  Diamond  Waveguide  (a)  Standard  single-pass  laser 
excitation  exhibiting  a  short  path  length  of  interaction  with  NV  centers  as  compared 
to  the  LTDW.  (b)  For  an  equivalent  excitation  power,  the  LTDW  increases  the  path 
length  of  the  excitation  beam  by  up  to  several  orders  of  magnitude,  (c)  The  colour 
CCD  image  of  the  LTDW,  excited  by  a  green  pump  laser  of  70  mW,  shows  bright 
fluorescence  without  any  spectral  filtering,  (d)  Simulation  of  the  beam  path  in  the 
LTDW  with  a  =  500  pm  input  facet,  3  mm  side  length,  100  pm  diameter  Gaussian 
profile  input  incident  at  11.6°,  and  an  absorption  constant  of  a  =  0.45  cm-1. 


Red  fluorescence  is  recorded  with  a  colour  CCD  camera  without  spectral  filtering. 
Figure  2-ld  shows  the  corresponding  numerical  simulation  of  a  Gaussian  pump  laser 
incident  on  an  LTDW  structure  with  the  same  dimensions.  For  clarity  in  visualizing 
light  propagation,  a  non-optimal  length  pump  path  of  only  6  cm  is  depicted. 

Simulations  of  path  length  versus  incident  angle  in  Fig.  2-2b  indicate  that  the 
LTDW  structure  enables  path  lengths  of  over  a  meter.  These  numerical  simulations 
(see  section  2.2.2)  account  for  beam  divergence  and  assume  losses  due  to  absorption 
only.  Scattering  losses  are  neglected  due  to  the  low  surface  roughness  of  the  polished 
sample  (see  section  2.5.2). 

In  this  ensemble  sensor,  the  available  pump  power  is  limited,  making  efficient 
conversion  of  pump  light  into  detected  ODMR  signal  essential.  We  quantify  this 
“pump-to-signal  photon  conversion”  as  1]odmr  =  (ftdet/(ftp,  where  (ft([et  ((ftp)  is  the  photon 
flux  detected  (in  the  pump  beam).  Figure  2-2a  depicts  (ftdet  vs.  (ftp  when  the  LTDW 
input  facet  is  pumped  with  a  532  nm  Gaussian  beam  (1/e2  diameter  =  300  pm).  Here, 
(ftp  is  obtained  from  the  power  of  the  pump  beam;  the  linear  relationship  indicates 
that  the  saturation  is  not  reached  even  at  relatively  high  pump  power  (~  2W).  The 
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Figure  2-2:  Optical  properties  of  the  LTDW  (a)  Pump  photon  flux  is  plotted 
against  detected  photon  flux  for  the  LTDW  (blue)  and  for  a  single  pass  through  the 
long  side  of  the  device  (red).  Error  bars  are  contained  within  the  data  points,  (b) 
Maximum  optical  path  length  (unconstrained  by  absorption)  is  plotted  as  a  function 
of  input  angle  into  the  LTDW  structure  with  a  150  pm  input  facet  and  standard 
beam  divergence  in  blue.  The  excitation  laser  intersects  the  input  facet  50  pm  from 
the  center,  resulting  in  asymmetry  in  the  plot.  For  comparison,  the  maximum  optical 
path  length  achievable  in  a  single-pass  configuration  through  the  top  of  the  sample 
is  also  plotted  in  black,  and  through  the  side  of  the  sample  is  plotted  in  red.  In 
samples  with  sufficiently  high  NV  density,  the  excitation  beam  would  be  absorbed 
before  achieving  the  maximum  path  length  allowed  by  the  LTDW  structure. 


slope  provides  a  conversion  efficiency  of  t]odmr  =  2.6  x  1CT3,  i.e.,  approximately 
one  red  photon  is  detected  for  every  385  green  pump  photons.  This  efficiency  is  two 
or  three  orders  of  magnitude  greater  than  when  the  diamond  is  pumped  in  a  single 
pass  configuration  along  the  long  or  short  axis  respectively.  From  our  measurements, 
the  fraction  of  pump  photons  converted  into  NV  fluorescence  (both  collected  and  not 
collected)  is  rjf  ~  0.055  (see  section  2.5.1). 


2.2.2  Matlab  and  Zemax  Simulation 

Several  key  factors  that  affect  the  path  of  the  input  beam  in  the  diamond  sample 
include  the  length  of  the  input  facet,  the  angle  the  beam  couples  into  the  sample,  the 
placement  of  the  beam  on  the  input  facet  as  it  enters  the  sample,  and  the  absorption 
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Figure  2-3:  (a)  A  MATLAB  simulation  example  of  a  pattern  where,  after  32  reflec¬ 
tions,  the  beam  exits  the  sample  at  an  angle  adequately  separated  from  the  input 
path  to  enable  transmission  measurements,  (b)  An  example  of  a  path  that  spatially 
covers  the  majority  of  the  diamond  sample,  (c)  A  three-dimensional  rendering  of 
the  LTDW  from  the  Zemax  simulation  featuring  a  diverging  Gaussian  beam.  Beam 
divergence  is  artificially  exaggerated  by  two  orders  of  magnitude  for  clarity. 


constant,  a,  of  the  sample.  For  the  MATLAB  simulations,  we  assume  that  polishing 
imperfections  at  the  edges  and  corners  of  the  sample  do  not  result  in  scattering 
loss.  Ray-tracing  models  are  expanded  to  illustrate  the  excitation  of  the  sample  by 
a  Gaussian  beam  shape  in  two  dimensions.  Zemax  simulations  consider  paraxial 
Gaussian  beams  and  assume  perfect  total  internal  reflection  and  losses  only  due  to 
material  absorption  in  a  three-dimensional  structure.  For  the  beam  diameter  in  our 
experimental  setup,  most  beam  paths  are  not  longer  than  the  Rayleigh  range.  To 
confirm  that  the  beam  continues  to  be  confined  even  when  the  Rayleigh  range  is 
much  shorter  and  to  confirm  that  interactions  with  the  top  and  bottom  face  of  the 
LTDW  would  not  result  in  beam  paths  exciting  the  device,  we  artificially  increased 
the  beam  divergence  by  two  orders  of  magnitude  to  exaggerate  the  possible  effect. 
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Even  at  this  extreme,  all  beam  paths  were  confined  until  interacting  with  the  input 
facet  at  an  angle  not  meeting  the  TIR  condition.  Paths  that  continue  after  partial 
reflections  were  not  considered. 

The  following  parameters  affect  path  length  in  the  LTDW: 

•  Input  Facet  Length:  For  maximizing  the  absorption  of  the  input  beam,  a  shorter 
input  facet  provides  fewer  opportunities  for  the  beam  to  out-couple  from  the 
structure.  A  longer  input  facet  is  preferable  for  transmission-type  experiments, 
as  it  allows  for  more  beam  path  configurations  where  the  beam  leaves  the  sam¬ 
ple  at  a  separate  angle  from  which  it  entered  -  practically  this  allows  room 
for  a  detector  to  be  placed  away  from  the  input  beam  path  for  transmission 
experiments. 

•  Input  Beam  Angle:  The  minimum  requirement  for  the  input  beam  angle  is  that 
it  must  satisfy  the  TIR  condition  on  both  its  first  and  second  reflection  in  the 
sample,  since  the  rectangular  cut  of  the  diamond  preserves  these  angles.  As  seen 
in  Figure  2- 2b,  incident  angles  that  are  more  than  ~  60°  from  incident  do  not 
satisfy  this  requirement  and  result  in  a  path  length  of  only  several  millimeters. 

•  Input  Beam  Placement:  If  the  input  beam  diameter  is  smaller  than  the  length 
of  the  input  facet,  placing  the  beam  away  from  the  center  of  the  input  facet 
breaks  the  symmetry  of  the  system  and  yields  a  wider  range  of  patterns. 

•  Rayleigh  Range:  The  high  index  of  diamond  (n  =  2.4)  aids  in  the  confinement  of 
the  beam,  however,  the  dispersion  of  the  input  beam  can  increase  the  interaction 
of  the  beam  with  edges  and  corners  that  are  more  likely  to  have  polishing 
imperfections,  resulting  in  scattering  losses.  A  more  tightly  focused  beam  can 
fit  into  a  smaller  input  facet,  but  will  have  a  shorter  Rayleigh  range. 

•  Absorption  Constant:  The  path  length  at  which  99%  of  the  beam  is  fully  ab¬ 
sorbed  by  the  sample  can  be  calculated  from  the  <x  value  for  a  given  wavelength 
for  the  sample. 
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2.3  Theory  and  Relevant  Calculations 


2.3.1  Estimation  of  number  of  NVs  addressed 


The  LTDW  sample  used  in  the  present  apparatus  has  a  volume  of  2.7  x  10“3  cm3.  A 
pure  carbon  lattice  of  this  volume  with  no  defects  corresponds  to  ~  5  x  1020  carbon 
atoms  (Eq.  2.1). 


2.7  x  10“3  cm3  x 


3.515  g 
cm3 


1  mol 
12.01  g 


6.02  x  1023  atoms 
1  mol 


4.76  x  102°  atoms  (2.1) 


The  absorption  constant  of  the  sample  (ct  ~  0.45  cm-1)  is  measured  by  single  pass 
beam  transmission  while  accounting  for  front  and  back  surface  reflections.  Assuming 
a  single  NV  absorption  cross  section  of  ~  2  x  10“ 17  cm2  [41],  this  yields  an  NV  density 
of  2.3  x  1016  cm  3  for  this  sample.  From  these  parameters,  the  NV  center  density  is 
estimated  at  0.1  ppm,  or  ~  1013  NV  centers  in  each  orientation.  This  is  the  maximum 
number  of  NV  centers  that  can  be  addressed.  Conservatively,  the  minimum  number 
of  NV  centers  addressed  is  estimated  to  be  two  orders  of  magnitude  smaller,  yielding 
more  than  five  orders  of  magnitude  improvement  over  a  single  NV  center  with  the 
same  excitation  and  collection  efficiency  parameters. 


2.3.2  Path  length  of  a  Gaussian  beam  in  the  LTDW 

The  path  length  of  the  excitation  beam  in  the  sample  is  calculated  geometrically  via 
ray-tracing  for  rays  of  incident  angle  6  with  respect  to  the  prism  facet,  in  the  plane 
of  the  diamond  waveguide.  For  a  given  incidence  location  on  the  input  facet,  this 
ray  gives  the  maximum  path  length  ( Lmax(6 ))  with  respect  to  the  input  angle.  We 
then  consider  a  Gaussian  beam  as  approximated  by  a  sum  of  rays.  In  particular,  we 
consider  a  Gaussian  beam  as  a  sum  of  rays  at  angles  6  with  weight  Fg(9,  0o,(j),  where 
sigma  is  the  FWHM  of  a  small  (NA<^t  1)  Gaussian  distribution  of  angles  around  6q. 
We  define  the  mean  path  length  of  this  Gaussian  beam  as  Lg{9q)  in  the  relation 
Pout,G  —  Pin,Ge~aLG^9°\  where  Pin  is  the  incident  beam  power,  Pout,G  is  the  beam 
power  exiting  the  diamond  after  some  pathlength  within  it,  and  a  is  the  absorption 
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constant  of  diamond  at  532  nm  for  a  given  NV  density.  We  assume  that  the  beam  is 
spatially  centered  at  the  middle  of  the  prism  facet.  Solving  for  Lq(9q )  gives 


Lg(9q)  =  —a  Mn 


=  —a  1  In 


(  Pout,G  \ 

V  Pin,G  ) 

(  f  FG(0,6o,a)e-aL™* Wd9\ 

V  f  FG(9,90,a)ffl  ) 


(2.2) 


If  the  NV  density  is  low  and  a  is  much  less  than  1/Lmax(9 )  —  i.e.,  the  losses  are 
dominated  by  scattering  rather  than  material  losses  —  we  approximate  e~aLmax ^  as 
1  —  aLmax(9)  and  ln(l  +  ()  as  (,  so  that  S2  becomes 


Lg(90) 


-i,  (}  ajFG(9,90,a)Lmax(9)d9\ 

a  nV  I  Fg(9,  90,a)d6  J 

(~ot  f  FG(9, 90,cr)Lmax(9)d9\ 

""  {  /FG(9,90,a)d9  J 

J  FG(9,90,a)Lmax(9)d9 
f  FG(9,90,a)d9 


(2.3) 


Equation  2.3  states  that  the  path  length  for  a  Gaussian  beam  in  the  LDTW  whose 
loss  is  dominated  by  scattering  (not  material  losses  inside  the  diamond)  is  a 
normalized,  Gaussian-weighted  average  of  the  path  lengths  calculated  in  the 
ray-tracing  approximation,  over  a  spread  of  angles  characterized  by  the  beam 
divergence  and  the  spatial  entry  of  the  beam.  In  other  words,  Lg(6q)  approximates 
the  longest  path  that  a  Gaussian  beam  of  angle  spread  o  and  central  angle  60 
propagates  in  the  LTDW. 


2.3.3  Effective  numerical  aperture  of  the  LTDW 

To  achieve  the  highest  number  of  reflections,  the  ideal  dimension  of  the  input  facet 
depends  on  the  coupling  method.  Since  beam  paths  that  refract  into  the  diamond 
through  the  input  facet  may  subsequently  couple  out  on  that  facet,  a  smaller  input 
facet  can  increase  the  confinement  of  light  inside  the  LTDW.  A  larger  input  facet 
allows  multiple  overlapping  TIR  patterns  to  occur  for  entrance  angles  which  reflect  off 
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the  facet  and  satisfy  the  TIR  condition  (6  >  9C  =  24.6°),  which  can  can  result  in  more 
uniform  excitation.  If  the  laser  is  focused  with  a  confocal  parameter  several  times  the 
diamond  width,  one  can  perform  transmission  measurements  using  a  photodiode  on 
the  beam  exiting  the  structure  through  the  coupling  facet.  The  minimum  requirement 
for  focusing  the  beam  is  that  the  beam  divergence  does  not  exceed  the  effective 
numerical  aperture  of  the  structure,  i.e.,  the  range  of  angles  that  satisfy  the  TIR 
condition  when  the  laser  is  incident  on  the  input  facet.  The  effective  NA  of  the 
symmetric  structure  used  here  is  0.84. 


2.3.4  Fraction  of  emission  from  an  emitter  deep  within  a  par¬ 
allelepiped 

A  disadvantage  to  the  high  refractive  index  of  diamond  is  that  a  large  portion  of  the 
NV  emission  is  trapped  in  the  diamond  and  won’t  exit  the  LTDW.  This  is  also  an 
issue  for  more  standard  confocal  collection  schemes.  For  this  calculation,  we  neglect 
re-absorption  and  assume  that  the  NV  emission  is  isotropic  -  this  is  not  the  case, 
however,  we  approximate  the  overlap  of  the  emission  from  all  four  orientations.  For  an 
emitter  deep  within  a  parallelepiped  of  refractive  index  n,  there  is  a  cone  defined  by  the 
critical  angle  where  the  emission  can  escape  each  facet  of  the  solid.  For  the  diamond- 
air  interface,  this  critical  angle  is:  9C  =  24.6.  The  ratio  of  the  solid  angle  where  the 
emission  can  escape  (defined  by  the  critical  angle  cone)  to  the  solid  angle  of  the  facet 
can  be  described  by:  27r(1V°h e°'1  =  3(1  —  cos(sin_1(^)))  «  3(1  —  \Jl  —  4j)  =  .272  or 
27.2%  of  the  emission  from  an  emitter  deep  within  a  parallelepiped  escapes  the  device 
and  can  be  detected  [44],  Taking  the  orientation  of  the  NV  dipole  into  account,  there 
may  be  an  advantage  to  side-collection;  we  point  to  the  supplementary  material  of 
ref.  [14],  which  states:  “The  photon  collection  efficiency  of  the  side-collection  method 
is  therefore  expected  to  be  in  the  range  of  29%  <  r]c  <  91%  if  many  diffuse  reflections 
between  the  edge  of  the  diamond  are  allowed”  and  goes  through  the  full  derivation 
taking  into  account  the  four  orientations  of  the  NV  centers. 


39 


2.4  Sample  preparation 


Sample  preparation  begins  with  a  (lOO)-oriented,  type  Ha  CVD  diamond  produced 
by  chemical  vapor  deposition  (Element  6)  with  dimensions  of  (3  x  3  x  0.3  mm3).  To 
increase  NV  ensemble  density,  the  diamond  was  electron- irradiated  at  4.5  MeV  with 
a  beam  current  of  20  mA  over  ~  0.15  m2  and  annealed  for  2  hours  at  850°C  yielding 
~  0.1  ppm  NV  centers.  Characterized  via  absorption  measurements  (a  ~  0.45  cm-1), 
this  concentration  has  adequate  brightness  to  illuminate  the  pattern  of  the  pump  laser 
in  the  sample  and  allow  for  almost  complete  absorption  of  the  pump  while  being  low 
enough  to  preserve  spin  coherence  properties.  At  low  pump  and  microwave  power 
levels,  the  ODMR  line  width  is  limited  by  the  coherence  of  the  diamond,  which 
is  ~  1  ps  as  characterized  by  free-induction  decay  measurements.  All  six  surfaces  of 
the  diamond  are  polished  to  a  mean  surface  roughness  of  ~  2  nm  RMS  (characterized 
via  atomic  force  microscopy,  see  section  2.5.2)  with  a  500  pm  input  facet  on  one 
corner  at  45°  creating  an  entrance  window  for  the  laser  (see  Fig.  2-la).  The  diamond 
sample  is  thermally  coupled  to  a  brass  sample  holder  with  indium,  which  may  cause 
additional  absorption  through  evanescent  coupling  of  the  confined  radiation.  The 
thermal  coupling  also  ensures  that  the  device  is  in  thermal  equilibrium  throughout 
the  experiments,  avoiding  broadening  in  the  ODMR  line  widths.  The  sample  holder 
is  temperature-stabilized  to  21.6  ±  0.1  C  using  a  water-chiller. 


2.5  Device  Characterization 

2.5.1  Conversion  Efficiency 

We  quantify  the  “pump-to-signal  photon  conversion”  as  tjodmr  =  ^det/^pi  where  (ftdet. 
(0p)  is  the  photon  flux  detected  (in  the  pump  beam).  Here,  (j)p  is  obtained  from 
the  power  of  the  pump  beam,  accounting  for  the  17%  intensity  reflection  due  to 
the  air/diamond  interface.  The  slope  from  figure  2-2a  in  the  main  text  provides 
the  conversion  efficiency  tjodmr  —  2.6  x  10~3,  i.e.,  approximately  one  red  photon 
is  detected  for  every  385  pump  photons.  This  efficiency  is  two  (three)  orders  of 
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magnitude  greater  than  when  the  diamond  is  pumped  in  a  single  pass  configuration 
along  the  long  (short)  axis.  We  can  also  deduce  from  our  measurements  the  fraction 
of  pump  photons  converted  into  NV  fluorescence,  i.e.,  rj /  =  where  the  total 

emitted  photon  flux  0/  is  obtained  from  the  fluorescence  detection  efficiency  i ]det  as 
0/  =  < Pdet/Vdet ■  We  estimate  rj(iet  ~  0.048  from  the  fraction  of  photons  emitted  into 
angles  that  are  not  confined  by  TIR  (27.2%  of  total  fluorescence)  and  the  fraction  of 
these  photons  that  are  imaged  onto  the  detector  (~  17%).  This  provides  r]f  ~  0.055. 
In  principle,  one  may  expect  ly  as  high  as  the  NV  internal  quantum  efficiency,  which 
has  been  approximated  to  be  near  0.7  [45]  for  NVs  in  bulk  diamond  due  to  the  charge 
state  transfer  to  NV°.  Because  we  were  not  able  to  accurately  measure  the  green  pump 
scattered  out  of  the  diamond,  it  was  not  possible  to  independently  estimate  the  NV 
quantum  efficiency  in  this  experiment.  However,  by  improving  diamond  polishing  to 
reduce  scattering  loss,  the  pump  beam  confinement  could  be  improved  to  increase  r// 
further.  The  detection  efficiency  could  also  be  increased  by  a  factor  of  2.5  by  directly 
positioning  detectors  on  the  four  diamond  facets,  as  was  done  in  Ref.  [14]. 

2.5.2  Surface  roughness  of  the  polished  diamond 

The  surface  roughness  of  the  polished  diamond  was  characterized  using  a  Digital 
Instruments  NanoScope  atomic  force  microscope  (AFM).  Figure  2-4  shows  a  3  pm  by 
3  pm  square,  representative  of  the  surface  features  present  across  the  sample.  The 
mean  surface  roughness  was  found  to  be  0.78  nm  and  the  rms  surface  roughness  was 
found  to  be  2.01  nm.  This  surface  roughness  is  skewed  by  sparse  pits  in  the  surface, 
which  could  be  ameliorated  with  improved  polishing  techniques. 

2.5.3  Absorption-limit  on  path  length 

In  the  case  of  the  diamond  shown  in  Fig.  2-1,  the  path  length  is  largely  limited  by  the 
absorption  of  the  beam.  Using  a  carefully  aligned  single-pass  beam,  we  calculate  the 
absorption  constant  of  this  sample  to  be  a  =  0.45  cm-1.  With  increased  NV  density 
we  would  expect  this  to  increase,  limiting  the  utility  of  the  LTDW  geometry  for  high 
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Figure  2-4:  Atomic  force  microscopy  characterization  of  surface  polishing, 
defect  density  samples. 

2.6  Outlook 

To  improve  on  the  LTDW  design  presented  here,  T1R  geometries  utilizing  quasi- 
chaotic  excitation  paths  could  also  be  used  for  NV  excitation  and  may  result  in  more 
uniform  coverage  of  the  bulk  sample  [46].  While  the  facets  of  the  LTDW  are  highly 
polished  to  a  surface  roughness  of  ~  2  nm  root-mean-square  (RMS)  (see  section  2.5.2), 
scattering  loss  occurs  at  the  edges  and  corners  of  the  structure.  Lower  surface  rough¬ 
ness  would  improve  the  TIR  and  reduce  scattering  losses.  Additionally,  increasing 
ODMR  contrast  also  would  result  in  improved  sensitivity.  Another  promising  geome¬ 
try  is  the  whispering  gallery  mode  resonator,  which  has  previously  been  demonstrated 
in  diamond  at  1R  frequencies  [47]. 
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Chapter  3 


Broadband  magnetometry  and 
temperature  sensing  using  a  LTDW 


Portions  of  this  work  appeared  in  print  here:  [39] 

3.1  Introduction 

The  LTDW’s  increased  path  length  and  the  consequently  large  ODMR  signal  photon 
flnx  allows  for  high  precision  measurements  of  quantities  that  affect  the  NV’s  mag¬ 
netic  sub-levels.  Figure  3-la  plots  ODMR  spectra  obtained  under  continuous-wave 
(CW)  microwave  held  excitation  when  a  3  mm  diameter  loop  of  200  pm- diameter 
wire  is  placed  5  mm  above  the  LTDW.  An  additional  2.5  mT  static  magnetic  held 
applied  along  the  (111)  diamond  crystal  axis  splits  the  degeneracy  of  the  NV  sub-level 
transitions  ms  —  0  — >■  ±1.  Of  the  four  NV  orientations  supported  in  the  diamond 
lattice,  the  (111)  sub-ensemble  therefore  shows  the  greatest  Zeeman  shift,  as  seen  in 
the  splitting  between  frequencies  cu+  and  o;_  (Fig.  3-la),  as  its  quantization  axis  is 
aligned  with  the  magnetic  held.  The  equal  magnetic  held  projection  onto  each  of 
the  three  remaining  NV  orientations  results  in  their  degenerate  transition  energies 
with  three  times  the  ODMR  contrast  compared  to  the  (111)  sub-ensemble.  Hyperhne 
coupling  of  the  NVs  to  the  14N  nuclear  spin  splits  each  of  the  electronic  sub-level  tran¬ 
sitions  into  Lorentzian  triplets  (Fig.  3c),  which  exhibit  a  full-width  at  half-maximum 
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linewidth  of  1.2  MHz  spaced  by  ~  2.1  MHz. 


Figure  3-1:  ODMR  spectra  and  level  structure  (a)  CW  electron  spin  resonance 
with  external  magnetic  field  aligned  along  a  single  (111)  crystal  axis.  The  inset  depicts 
the  diamond  lattice,  where  the  blue  circle  represents  the  nitrogen  atom  and  the  green 
circles  show  the  four  possible  orientations  for  the  vacancy  in  the  tetrahedral  crystal 
lattice,  constituting  the  four  sub-ensembles  of  NVs.  (b)  Lock-in  output  corresponding 
to  signal  in  (a).  The  scale  factor  (V/ Hz)  is  provided  by  a  linear  fit  around  each  of  the 
two  red  points.  Temperature  and  magnetic  held  shifts  are  measured  independently 
by  tracking  both  ms  =  ±1.  (c)  Energy-level  diagram  of  diamond  NV  center  showing 
radiative  (solid  lines)  and  non-radiative  (dotted  lines)  transitions. 


3.2  Sensitivity  measurements 

3.2.1  Noise  floor 

Frequency  analysis  of  the  noise  in  the  LTDW  system  (Fig.  3-2a)  indicates  that 
at  1  Hz,  one  side  of  the  (111)  orientation  sub-ensemble  resonance  is  sensitive  to 
less  than  ~  1  iiT/a/Hz  (Fig.  3-2a).  However,  as  both  sides  of  the  resonance  are 
not  being  monitored,  the  bulk  of  the  noise  in  this  measurement  can  be  attributed 
to  thermal  fluctuations  in  the  environment.  We  expect  the  next  most  prominent 
noise  contribution  to  be  electronic  noise  from  the  system,  followed  by  environmental 
magnetic  held  noise,  as  this  setup  was  not  protected  by  magnetic  shielding. 
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3.2.2  Separation  of  magnetic  field  and  temperature  data 


We  measure  magnetic  field  and  temperature  shifts  by  monitoring  the  fluorescence  of 
the  LTDW  at  the  ms  =  0—^—1  and  the  ms  =  0  — )•  +1  transitions.  The  microwave 
excitation  frequency,  cuodmr,  is  alternately  tuned  to  the  frequencies  uj+  and  (u_,  where 
the  derivative  of  the  ODMR  signal  with  respect  to  frequency  is  greatest.  Around  these 
frequencies,  we  modulate  cuodmr  at  1.5  kHz  with  a  modulation  depth  of  1  MHz  and 
use  a  lock-in  amplifier  to  monitor  the  NV  fluorescence.  Figure  3-lb  depicts  the  lock-in 
signal  for  a  time  constant  of  10  ms. 


Frequency  (Hz) 


Figure  3-2:  Sensitivity  measurements  (a)  Noise  density  of  the  in-phase  lock-in 
amplifier  signal  at  the  cu+  resonance  (see  inset)  converted  to  units  of  magnetic  held 
for  a  time  constant  of  100  ms;  low  frequency  noise  is  attributed  to  thermal  effects. 
(b)ESeparation  of  magnetic  held  and  temperature  effects.  Green  and  blue  lines  show 
the  sum  and  difference  of  the  intensity  signals  from  lower  and  higher  frequency  res¬ 
onances  of  the  (111)  orientation  with  a  time  constant  of  10  ms,  decoupling  the  tem¬ 
perature  and  magnetic  held  drifts  over  100  s,  respectively.  A  20-s-period  square  wave 
of  amplitude  190  nT  is  applied  with  external  coils  and  is  clearly  visible  in  the  blue 
trace.  The  noise  on  a  single  10  s  step  (0.93  nT)  extrapolates  to  ~  290  pT  Hz-1/2 
when  accounting  for  the  10  Hz  equivalent  noise  bandwidth  of  the  lock-in  amplifier. 
The  green  trace  shows  drift  attributed  to  environmental  effects;  aliasing  effects  have 
been  removed  (less  than  200  ms). 


From  the  ground-state,  axial,  hyperhne  Hamiltonian  (see  section  1.1.2),  it  follows 
that  the  ms  =  0  — >■  ±1  resonance  frequencies  are  given  by  cj±  ~  ( Dgs  +  /3tAT)  ±~/Bz, 
where  (3t  ~  —74  kHz/K  is  the  phenomenological  temperature  dependence  of  the 
NV  center’s  zero  held  splitting  (ZFS)  at  room  temperature  [33]  and  the  Zeeman 
shift  in  response  to  an  axially  applied  magnetic  held  in  the  (111)  crystal  direction  is 
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7  =  \^Bds/h  ~  28  GHz/T,  where  \1b  =  (eft)/(2me)  is  the  Bohr  magneton,  and  gs  ~  2 
is  the  Lande  g- factor.  By  sequentially  monitoring  both  transitions,  oj±,  we  obtain 
the  mean  value  of  uj±,  which  yields  the  temperature  shift  AT,  and  the  difference, 
uj+  —  ta_,  which  provides  the  magnetic  held  shift  [48]. 

Figure  3-2b  plots  independent  measurements  of  temperature  and  magnetic  held 
over  a  time  period  of  100  seconds.  Each  point  represents  an  integration  time  of  10  ms. 
We  resolve  sub-nT-level  noise  on  a  190  nT  magnetic  held  applied  as  a  0.05  Hz  square- 
wave  and  expect  a  decrease  in  noise  for  integration  times  of  one  second.  The  noise  on  a 
single  10  s  step  (0.93  nT)  extrapolates  to  ~  290  pT  _1/2  when  accounting  for  the  10  Hz 
equivalent  noise  bandwidth  of  the  lock-in  ampliher.  We  separately  resolve  a  gradual 
temperature  shift  during  this  measurement  (green  curve),  within  the  expected  range 
of  temperature  stabilization.  Note  that  for  this  method  of  temperature  and  magnetic 
held  shift  separation,  signal  errors  add  in  quadrature. 

3.2.3  Simultaneous  measurement  of  Zeeman-split  resonances 

Unlike  the  work  described  above,  which  utilizes  sequential  measurement  of  Zeeman- 
split  resonances,  a  frequency-locking  technique  was  also  developed  to  monitor  both 
resonances  simultaneously.  A  closed-loop  system  that  locks  to  the  center  frequency 
of  an  NV  resonance  allows  for  the  steady-state  operation  of  an  NV  magnetometer  to 
be  isolated  from  the  phenomenological  variables  that  determine  the  scale  factor  of 
previous  lock- in-based  approaches  [49,  50,  51].  Furthermore,  by  virtue  of  locking  to 
the  NV  resonance  frequency,  measurement  is  always  performed  in  the  approximately 
linear  regime  of  the  lock-in  signal,  and  the  full  dynamic  range  of  the  NV  center  can 
again  be  accessed.  The  design  has  similar  benefits  and  drawbacks  as  an  atomic  Mz 
magnetometer  [52] ;  in  particular,  the  bandwidth  of  an  Mz  conhguration  is  limited  by 
the  modulation  frequency  which  must  be  less  than  1  /T\  where  T\  is  the  longitudinal 
spin  relaxation  rate.  The  NV  T\  under  optical  excitation  limits  the  modulation  fre¬ 
quency  to  <  20  kHz  [49,  51].  Via  frequency-generation  and  lock-in  implementation 
on  a  custom  built  FPGA  board,  all  four  sub-ensembles  can  be  addressed  for  vector 
measurements  of  magnetic  fields. 
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The  combination  of  atom-like  long-term  stability  and  vector  capability  in  a  com¬ 
pact,  solid-state  package  makes  NV-based  magnetometers  a  prime  candidate  for  use  in 
a  wide  variety  of  field  applications,  from  all  magnetic  navigation  [53,  54]  to  prospect¬ 
ing  for  oil  wells  [55]  to  magnetic  mapping  of  archaeological  sites  [56]  to  locating 
unexploded  ordinance  [57,  58,  59].  For  practical  operation  in  an  unshielded  field  en¬ 
vironment,  however,  a  vector  magnetic  sensor  also  requires  high  dynamic  range  to  be 
robust  against  rotation  within  the  Earth’s  magnetic  field,  which  may  sweep  measured 
magnetic  field  projections  over  a  range  up  to  ~  130  pT.  Using  a  frequency  locking, 
multi-axis  technique,  scale-factor  free,  temperature-independent,  vector  magnetom- 
etry  is  demonstrated  over  a  high  dynamic  range.  This  is  a  vital  step  forwards  in 
taking  the  NV  magnetometer  from  a  laboratory  demonstration  to  a  functional  device 
for  detecting  fields  in  an  unshielded  environment. 


3.3  Experimental  details 

3.3.1  Apparatus 

A  spatially  filtered,  532  nm  single-mode  longitudinal  excitation  beam  from  a  Verdi 
V5  or  a  homebuilt  50  mW  direct  green  laser  diode  (OSRAM)  is  focused  into  the  input 
facet  of  the  diamond  sample.  Both  frequency  and  spatial  multi-mode  excitation  is 
feasible  as  long  as  the  beam  shape  can  be  focused  to  the  size  of  the  input  facet  i.e., 
below  300  pm  diameter.  Three  sets  of  perpendicular  Helmholtz  coils,  with  diameters 
ranging  from  8  to  12  inches,  provide  a  uniform  static  or  slowly  varying  magnetic 
held.  Noise  floor  measurements  were  taken  with  a  Halbach  ring  to  provide  a  static 
magnetic  held.  Using  ODMR  spectra  recorded  with  the  application  of  three  different 
control  helds,  the  rotation  matrix  from  the  cartesian  frame  of  the  coils  to  the  NV 
axis  is  determined.  An  impedance-matched  loop  antenna  located  ~  5  mm  above  the 
sample  delivers  microwave  excitation,  with  40  dBm  on  the  transmission  line.  The 
microwave  excitation  is  frequency  modulated  at  1.5  kHz  with  a  modulation  depth  of 
1  MHz.  Fluorescence  emission  is  collected  with  an  aspheric  condenser  lens  (NA=0.79), 
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followed  by  a  662  —  800  nm  bandpass  filter  to  minimize  NV°  contribution  to  the 
fluorescence.  A  second  condenser  lens  images  the  fluorescence  onto  a  balanced  large- 
area  photoreciever  (New  Focus  2307),  which  minimizes  the  effect  of  low  frequency 
laser  noise  by  balancing  green  input  laser  power.  The  system  collection  efficiency 
is  ~  17%  of  the  total  fluorescence  signal,  which  is  demodulated  using  an  SRS-850 
lock-in  amplifier. 

3.3.2  Advantages  of  using  lock-in  amplification 

The  use  of  lock-in  amplifications  has  several  advantages  over  direct  readout  of  the 
ODMR  curve.  First,  by  modulating  and  demodulating  the  signal,  the  bulk  of  low 
frequency  electronics  noise  can  be  removed  from  the  signal.  This  is  analogous  to  the 
advantage  gained  by  using  dynamical  decoupling  techniques  to  filter  measurements  of 
specific  target  frequencies,  however  these  pulsed  techniques  are  also  able  to  provide 
further  improvement  in  sensitivity  by  taking  advantage  of  the  NV  coherence.  Sec¬ 
ondly,  the  lock-in  amplifier  produces  a  derivative  signal  which  can  be  approximated  to 
be  linear  as  it  crosses  zero.  This  vastly  simplifies  the  curve  fitting  necessary  for  mea¬ 
surements,  as  the  slope  of  the  line  can  be  measured  initially  and  checked  periodically 
(hourly  or  daily)  instead  of  fitting  Lorentzians  to  every  acquired  curve.  To  perform 
real-time  sensing,  extensive  curve-fitting  steps  need  to  be  minimized.  For  simulta¬ 
neous  measurement  techniques  utilizing  frequency-locking,  curve-fitting  is  eliminated 
altogether. 

3.3.3  Static  magnetic  field  alignment 

The  magnetic  fields  were  aligned  to  the  NV  axis  using  the  iterative  application  of 
incremental  test  fields  using  low-noise  current  sources  and  a  three-axis  Helmholtz  coil. 
Mat  lab  was  used  to  generate  a  rotation  matrix  from  the  lab  frame  to  the  NV  frame. 
For  the  work  presented  in  this  chapter,  we  confirmed  the  alignment  by  minimizing 
the  linewidth  of  the  other  three  sub-ensembles  -  when  the  total  magnetic  held  is 
aligned  with  one  orientation,  the  other  three  orientations  experience  a  degenerate 
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field  projection. 


3.4  Conclusions  and  Outlook 

3.4.1  Comparison  with  other  magnetometers 

Figure  3-3  shows  the  performance  of  the  LTDW  sensor  compared  to  recently  demon¬ 
strated  NV  magnetometers  and  other  magnetometer  platforms.  The  system  noise 
limit  from  Figure  3-2a  is  superimposed  over  the  frequency  range.  Also  indicated 
are  the  sensitivity  requirements  for  three  key  magnetometer  applications:  magneto- 
cardiography,  magnetoencelephography,  and  the  monitoring  of  geomagnetics.  The 
unique  position  of  the  LTDW  in  the  low-frequency  limit  makes  it  suited  for  a  variety 
of  magnetic  field  sensing  applications  including  medical  monitoring,  object  detection, 
and  the  study  of  geomagnetics  which  occur  in  the  10“4  Hz  to  10  Hz  range.  While 
SQUID  and  atom-based  magnetometry  achieve  better  sensitivity  at  low  frequencies, 
these  technologies  require  cryogenics  and  vacuum  systems  that  inhibit  portability  and 
versatility. 

3.4.2  Fundamental  spin  projection  limit 

For  high-frequency  magnetic  field  measurements,  Echo-type  magnetometry  could  be 
employed  with  dynamic  decoupling  pulses  to  extend  the  NV  electron  spin  coherence 
time  [8].  In  this  regime,  the  sensitivity  could  approach  the  fundamental  spin  projec¬ 
tion  limit, 

S  =  SBminVf  ~  (3.1) 

7  v  Nt 

where  5Bmm  is  the  minimum  detectable  magnetic  field,  T  is  the  measurement  time, 
N  =  1013  is  the  estimated  total  number  of  NV  centers  in  our  sample  (see  section 
2.3.1),  and  r  ~  1  ms  is  the  electron  spin  dephasing  time  measured  with  dynamical 
decoupling  in  comparable  samples  and  expected  to  be  feasible  for  NV  densities  on 
the  order  of  1  ppm  [8].  Neglecting  contrast  and  collection  efficiency,  these  parameters 
predict  a  spin-projection  limit  of  S'  =  0.36  fT/\/Hz. 
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Figure  3-3:  Comparison  of  the  LTDW  (dark  blue  overlay)  with  other  mag- 
netometry  technologies  Particularly  with  regard  to  target  applications  such  as  ge¬ 
omagnetics  [1]  shown  in  green  (corresponding  to  the  right  axis),  the  LTDW  achieves 
improved  sensitivity  with  respect  to  other  diamond-based  sensors,  shown  in  blue  (cor¬ 
responding  to  the  left  axis),  in  the  important  low-frequency  band  [2],  The  LTDW’s 
performance  is  competitive  with  other  magnetometer  platforms  [3],  shown  in  light 
blue. 


3.4.3  Next  steps  and  future  applications 

We  have  introduced  a  dual  magnetic  field  and  temperature  sensor  based  on  a  light¬ 
trapping  diamond  waveguide  geometry  that  enables  nearly  complete  green  pump  ab¬ 
sorption  and  efficient  spin-dependent  ffuorescence  collection.  With  more  than  two 
orders  of  magnitude  improvement  in  pump  power  conversion  over  previous  single¬ 
pass  schemes,  the  LTDW  achieves  high  signal  to  noise  sensing  and  sub-nT/VHz 
sensitivity  in  the  low-frequency  regime,  which  is  not  currently  accessible  by  dynami¬ 
cal  decoupling  techniques.  The  low  power  requirement  for  the  LTDW  allows  it  to  be 
implemented  with  a  compact,  inexpensive  laser  diode  excitation  source.  We  measure 
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a  noise  floor  of  less  than  300  pT / v^Hz  in  the  sub-Hz  regime.  In  this  low-frequency 
domain,  such  high-sensitivity  magnetometers  promise  new  applications  for  solid-state 
sensors  in  a  range  of  fields  including  biomedical  sensing  and  object  detection.  The 
light  trapping  technique  could  be  extended  to  transmission  measurements  and  di¬ 
rect  IR  absorption  techniques  [15].  The  device  is  not  limited  to  sensing;  it  could  be 
used  for  any  application  requiring  a  long  optical  path  length  within  the  sample  such 
as  electromagnetically-induced  transparency-based  experiments  or  optical  quantum 
memories  [60].  The  LTDW  provides  a  compact,  portable  precision  sensor  platform 
for  measuring  magnetic  fields,  temperature,  pressure,  rotation,  or  time.  Via  multi¬ 
plexing,  combinations  of  these  sensing  applications  could  be  performed  in  a  single 
compact  device. 
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Chapter  4 


Diamond  NV  electronic  and  nuclear 
spin-state  anti-crossings  under  weak 
transverse  fields 


This  chapter  is  adapted  from  [61],  which  was  written  in  close  collaboration  with 
Edward  H.  Chen. 


4.1  Introduction 

We  report  on  detailed  studies  of  electronic  and  nuclear  spin  states  in  the  diamond 
nitrogen  vacancy  (NV)  center  under  weak  transverse  magnetic  fields.  We  numerically 
predict  and  experimentally  verify  a  previously  unobserved  NV  hyperfine  level  anti¬ 
crossing  (LAC)  occurring  at  bias  fields  of  tens  of  gauss  -  two  orders  of  magnitude 
lower  than  previously  reported  LACs  at  ~  500  G  and  ~  1000  G  axial  magnetic  fields. 
We  then  discuss  how  the  NV  ground  state  Hamiltonian  can  be  manipulated  in  this 
regime  to  tailor  the  NV’s  sensitivity  to  environmental  factors  and  to  polarize  the 
nuclear  spin  state. 


55 


4.2  Energy  levels  of  the  diamond  NV  in  the  weak 
transverse  magnetic  field  regime 

4.2.1  Previous  work  in  the  transverse  magnetic  field  regime 
and  with  level  anti-crossings 

NV  defect  centers  in  diamond  are  optically  polarizable  quantum  systems  with  spin- 
dependent  fluorescence.  Using  electron  spin  resonance  (ESR)  under  ambient  con¬ 
ditions,  sensitivity  to  electric  fields  [16,  62,  63,  64],  transverse  and  axial  magnetic 
fields  [7,  8,  9,  10,  11,  65],  temperature  [18,  19,  20,  35],  strain  [66],  and  pressure  [67] 
have  been  observed  via  resonance  frequency  shifts  of  the  NV  ground-state  manifold. 
Non-separable  sensitivity  to  multiple  environmental  factors  is  problematic  for  to  using 
the  NV  as  a  sensor.  However,  the  Hamiltonian  governing  the  measurable  frequency 
shifts  can  be  tailored  to  enhance  (or  to  suppress)  sensitivity  to  different  physical  phe¬ 
nomena.  A  magnetic  bias  field  applied  parallel  or  perpendicular  (B ||  or  Bj_  as  shown 
in  Fig.  4- la)  to  the  NV’s  axis  in  the  diamond  crystal  lattice  energetically  separates 
the  spin  states  and  increases  sensitivity  to  magnetic  or  electric  fields,  respectively. 

Here,  we  investigate  an  unexplored  weak-field  regime  in  which  electronic  spin 
ground  state  energy  level  splittings  are  on  par  with  the  Zeeman  shift  induced  by 
an  applied  magnetic  field.  We  account  for  both  the  electron  and  the  nuclear  spin 
of  the  NV,  which  reveals  complex  dynamics  of  nuclear  spin  state  degeneracy  and 
previously  unobserved  hyperfine  level  anti-crossings.  These  features  occur  at  low 
magnetic  field  (B±  <  40  G)  as  compared  to  the  B\\  ~  500  G  and  B\\  ~  1000  G 
excited  and  ground  state  crossings  [68,  69,  70,  71],  which  have  been  used  for  nuclear 
spin  polarization,  providing  increased  sensitivity  to  resonance  shifts  through  narrower 
effective  linewidth  and  increased  contrast  [72,  36,  73,  74],  We  find  excellent  agreement 
between  experiment  and  theory  and  discuss  the  utility  of  the  nuclear  spin  degeneracy 
regime  toward  NV  sensing  applications  and  solid-state  atomic  memories  based  on 
nuclear  spin  polarization.  While  the  results  described  here  are  specific  to  the  NV, 
similar  anti-crossings  are  expected  in  any  spin-1  (or  higher)  defect  center  that  shows 
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hyperfine  level  splitting  on  the  same  order  of  magnitude  as  double-electron  spin  flip 
anti-crossings. 


Figure  4-1:  (a)  Axial  and  transverse  magnetic  field  directions  with  respect  to  NV 
axis.  As  B±  is  increased  from  35  G  (b)  to  45  G  (c),  the  mj  =  0  energy  levels  are 
separated  from  the  m;  =  ±1  energy  levels,  removing  the  degeneracy  which  results 
in  an  electron- nuclear  spin  flip  anti- crossing,  which  is  detailed  in  the  inset  of  (b). 
Non-dressed  states  are  indicated  with  colored  lines  and  dressed  states  are  indicated 
with  black  lines. 


4.2.2  Detailed  Hamiltonian 

The  NV  is  a  two-site  defect  with  a  spin-1  electronic  ground  state,  which  is  magnetically 
coupled  to  nearby  nuclear  spins.  For  a  single  NV  orientation,  energy  level  shifts  are 
described  by  the  following  spin  Hamiltonian  of  the  ground  triplet  state  in  the  presence 
of  magnetic,  electric,  and  strain  fields  [27],  taking  into  account  the  zero-held  splitting, 
nuclear  and  electronic  Zeeman  shifts,  Stark  shifts,  hyperfine  splitting,  and  nuclear 
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quadrupole  effects: 


Ugs  —{hDgs  +  d\\Uz)(S2z  -  i S(S  +  1)) 

-  d±[Ux(S2x  -  S2)  +  ILy(SxSy  +  SySX)} 

(4-1) 

I^BQe^S  •  -£?)  fingn(B  •  /) 

+  A\\SZIZ  +  A±(SXIX  +  Syly)  +  Pgs{lz - ^“) 

where  hDgs  is  the  NV  ground  state  zero  held  splitting  energy  (which  is  temperature- 
dependent  [35]),  d_ l  and  d\\  are  the  components  of  the  ground  state  electric  dipole 
moment,  the  total  effective  electric  held  H  —  E  +  a  encompasses  both  static  electric 
helds  ( E )  and  strain  (a),  ge  and  gn  are  the  electric  and  nuclear  Lande  g-factors,  Hb 
and  /j,n  are  the  Bohr  and  nuclear  magneton  constants,  B  is  the  applied  magnetic  held, 
A\\  and  Aj_  describe  the  axial  and  transverse  magnetic  hyperhne  interactions  with  the 
14N  nucleus,  Pgs  is  the  nuclear  electric  quadrupole  parameter,  S  is  the  electron  spin 
operator,  and  I  is  the  spin  operator  of  the  14N  nucleus. 

4.2.3  Double  electron  spin  flip  anti-crossing 

Figure  4-1  shows  the  numerically  calculated  electronic  ground  state  triplet  and  14N 
hyperhne  energy  levels,  \ms,mi),  where  rns  is  the  electronic  spin  state  and  mi  is 
the  nuclear  spin  state,  as  a  function  of  axial  magnetic  held  for  a  hxed  transverse 
magnetic  held.  State-mixing-induced  anti-crossings  are  seen  when  the  dressed  states 
(black  lines)  do  not  follow  high-axial- held  eigenstates  (solid  and  dashed  colored  lines). 

With  the  addition  of  a  weak  (<  3  G)  axial  magnetic  held,  we  see  mixing  and 
crossing  of  the  energy  levels,  resulting  from  off-diagonal  terms  in  1-Lgs  in  the 
basis.  The  double-electron  spin  hip  anti-crossings  occur  as  the  non-dressed  ms  =  1 
states  cross  ms  =  —1  states  with  identical  nuclear  spin,  for  example,  as  1 1, 0)  crosses 
|  —  1,  0) .  The  transverse  magnetic  held  leads  to  second-order  mixing  of  the  ms  = 
±1  states,  therefore  the  coupling  strength  Eg  of  these  electronic  Zeeman  interaction 
driven  level  anti-crossings  scales  quadratically  with  the  applied  transverse  magnetic 


field  as  described  in  [27]. 


4.2.4  Electron-nuclear  spin  flip  anti-crossing 

For  transverse  magnetic  fields  under  ~  40G,  as  shown  in  Fig.  4-lb,  an  additional  level 
anti-crossing  arises.  As  this  phenomenon  has  not  to  our  knowledge  been  previously 
described,  we  investigate  it  here  in  detail.  The  dressed  rn j  =  0  levels  (black  lines 
roughly  following  the  red  dashed  and  solid  lines)  cross  the  dressed  rn  j  =  ±1  levels 
(black  lines  roughly  following  the  blue  and  yellow  dashed  and  solid  lines),  resulting 
in  further  state  mixing  and  an  electron-nuclear  spin  flip  anti-crossing.  Solving  for  the 
eigenstates  of  7-Lgs  indicates  that  transverse  magnetic  fields  induce  state  mixing  in 
the  low  axial  magnetic  field  regime  by  bringing  the  nuclear  spin  states  into  resonance 
via  the  electronic  Zeeman  interaction  \fj,Bge(S  ■  B )]  in  the  presence  of  the  transverse 
hyperfine  interaction  [A±(SXIX  +  SyIy )].  This  state  mixing  results  in  dressed  states 
and  affects  which  transitions  are  allowed  by  optical  transition  selection  rules.  As  only 
electron  spin  flips  ( Sx  and  Sy  operators)  are  considered  in  our  calculation,  the  nuclear 
spin  flip  results  from  driving  the  electron  spin  flip.  This  is  a  spin-conserving  process. 
For  transverse  magnetic  fields  larger  than  ~  40  G,  the  energy  separation  added  by  the 
transverse  magnetic  field  removes  the  rn  j  =  0,  mj  =  ±1  level  degeneracy;  therefore, 
this  additional  anti-crossing  is  no  longer  present,  leaving  only  the  double-electron  spin 
flip  anti- crossing,  which  is  not  mediated  by  the  hyperfine  interaction. 

Figure  4-2  depicts  the  double-electron  spin  flip  and  electron- nuclear  spin  flip  anti¬ 
crossing  coupling  strengths  ( Eg ,  Ng)  and  the  axial  magnetic  field  positions  ( Ep ,  Np) 
of  these  anti-crossings  as  a  function  of  transverse  magnetic  field.  As  predicted  in  Ref. 
[27],  the  coupling  strength  of  the  double-electron  spin  flip  anti-crossing  is  proportional 
to  the  applied  transverse  magnetic  field  squared:  Eg  ~  (/isge.B_|_)2 / (Dgs  +  d||IIz)  in 
the  weak  electric  field,  weak  axial  magnetic  field  regime.  These  electron-spin  flip  anti¬ 
crossings  are  centered  around  Ep  =  0,  ±A\\h/ gefiB-  While  Eg  continues  to  increase  as 
B  l  increases,  Ng  vanishes  at  B±  ~  40  G  when  the  energy  levels  no  longer  overlap. 
The  axial  magnetic  field  position  about  which  this  anti-crossing  is  centered  (Np)  goes 
to  zero  with  increasing  B±  (see  App.  B). 
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B±  (G)  B_ l  (G) 

Figure  4-2:  (a)  Anti-crossing  coupling  strength  with  increasing  B±.  (b)  Magnitude  of 
small  axial  B  fields  at  which  the  centers  of  these  anti-crossings  are  observed.  In  both 
plots,  the  double-electron  spin  flip  anti-crossing  is  marked  E  and  the  electron-nuclear 
spin  flip  anti-crossing  is  marked  N.  At  fields  greater  than  B±  ~  40  G,  anti-crossing 
E  continues  to  increase  quadratically  and  anti-crossing  N  goes  to  zero  as  the  energy 
levels  no  longer  overlap. 


Using  second  order  degenerate  perturbation  theory,  we  calculate  the  change  in 
energy  as  a  function  of  applied  magnetic  and  electric/strain  fields.  By  solving  for  the 
axial  magnetic  held  (U||)  at  which  the  intersection  of  the  | mj  =  0)  and  the  | mj  =  1) 
energy  levels  occur,  we  obtain  a  closed-form  solution  for  the  location  of  the  electron- 
nuclear  spin  hip  level  anti-crossing.  This  is  in  good  agreement  with  our  numerical 
results  as  seen  in  Fig.  4-2. 


Np  = 


-a,d„(  +  \p,s\  J((d^  +  bi 

2  DU 


(4.2) 


where  £  =  Ay  —  P^s,  A\\  is  the  axial  hyperhne  interaction,  Dgs  is  the  zero  held  split¬ 
ting,  Pgs  is  the  quadrupole  energy,  and  Bj_  transverse  component  of  the  magnetic 
held  in  the  NV  frame  of  reference.  Figure  4-5  depicts  this  transverse  held  regime. 
Numerical  solutions  to  accounting  for  selection  rules  are  shown  in  conjunction 
with  experimental  data. 
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4.3  Experimental  apparatus 
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Figure  4-3:  (a)  Experimental  setup,  (b)  The  four  orientation  sub-ensembles  of  NV 
centers  are  shown  with  an  applied  magnetic  field,  (c)  An  anti-crossing  strength  of 
~  2.6  MHz  is  seen  in  sub-ensemble  1  when  the  perpendicular  component  of  the  B 
held  is  equal  to  30  G  and  the  axial  component  of  the  B  held  (iJqi)  is  varied 

around  zero  held.  Inset  shows  a  high  resolution  spectra  of  resonance  1  and  V . 


We  experimentally  validate  this  anti-crossing  by  performing  ESR  measurements 
on  an  NV  ensemble.  The  epitaxially-grown  diamond  is  polished  to  trap  the  green 
pump  light  and  guide  the  red  spin-dependent  huorescence  to  a  photodetector  in  a 
light-trapping  diamond  waveguide  geometry  [39].  We  estimate  the  intrinsic  strain 
in  this  sample  to  be  ~  3  x  10~5  from  the  ~  600  kHz  strain  splitting  of  the  central 
hyperhne  resonance  at  zero  applied  held,  given  values  of  d±  from  the  literature  [17]. 
We  expect  this  to  be  the  average  strain  experienced  across  the  sample.  The  NV 
has  four  orientations  in  the  diamond  lattice,  labeled  as  k  =  1,2,  3, 4  as  seen  in  Fig. 
4-3b.  We  lift  the  ms  =  ±1  degeneracy  and  the  orientation  degeneracy  by  applying 
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a  static  magnetic  field  at  an  angle  resulting  in  non-equivalent  projections  onto  each 
orientation.  A  60  G  Halbach  array  (magnitude  uniformity  >  99%  over  a  2  cm3 
volume  at  the  center)  is  positioned  around  the  sample  (Fig.  4-3a)  to  provide  this 
static  magnetic  held.  In  addition,  up  to  ±30  G  can  be  applied  in  x,  y,  or  5  (lab 
frame)  using  Helmholtz  coils.  The  total  magnetic  held  is  aligned  to  be  perpendicular 
to  the  k  —  1  orientation.  A  weak  axial  magnetic  held  sweep  (-E%||)  in  addition  to  this 
transverse  held  confirms  the  anti-crossings  in  the  ESR  spectra  (Fig.  4-3c).  Each  NV 
orientation  produces  two  Zeeman-split  34N  nuclear  hyperhne  triplets.  These  pairs 
are  located  symmetrically  around  c uc  =  Dgs  ±  Bj_.  In  levels  1  and  1',  for  which 
the  magnetic  held  is  perpendicular,  we  observe  both  predicted  anti-crossings  around 
Bi  ||  =  Np.  This  is  highlighted  in  the  inset  to  Fig.  4-3c. 

4.4  Experimental  verification 

4.4.1  Double  electron  spin  flip  anti-crossing 

While  Fig.  4-1  gives  insight  to  the  origin  of  the  level  anti-crossings  in  the  energy 
level  diagram,  Fig.  4-4  shows  the  corresponding  energy  level  transitions  at  B±  = 
35  G  and  45  G.  Numerical  solutions  to  T-Lgs  accounting  for  Sx  and  Sy  operators 
in  the  selection  rules  are  plotted  in  conjunction  with  the  experimental  data.  For 
both  values  of  transverse  held,  we  clearly  see  double-electron  spin  hip  anti-crossings 
centered  at  B\\  =  0,  ±A\\h/gefj,B  G  (Fig.  4-2b)  with  anti-crossing  coupling  strengths 
corresponding  to  those  expected  from  Fig.  4-2a.  As  expected,  we  also  see  the  center 
frequency  between  the  split  resonances  increase  with  increased  applied  transverse 
magnetic  held. 

4.4.2  Electron-nuclear  spin  flip  anti-crossing 

Figure  4-5  concentrates  on  the  region  of  interest  at  Bit±  =  30  G  around  which  the 
electron- nuclear  spin  hip  anti-crossing  is  maximized.  We  observe  the  electron-nuclear 
spin  hip  anti-crossing  at  ~  ±0.35  G  applied  axial  magnetic  held  as  predicted  by 
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Figure  4-4:  Plots  of  transition  energies  as  a  function  of  axial  magnetic  field  under 
fixed  transverse  magnetic  field.  Transitions  for  B  j_  =  35  G  and  B±  =  45  G  are  shown, 
which  correspond  to  the  energy  levels  in  Fig.  4-1.  Greyscale  indicates  normalized 
fluorescence.  Double-electron  spin  flip  anti-crossings  of  coupling  strength  Eg  are  seen 
at  Ep  =  0,  EA\\h/ geHB  —  ±0.75  G.  ESR  in  diamond  with  a  natural  abundance 
of  carbon  isotopes  has  additional  resonances  due  to  hyperfine  interactions  with  the 
~  1%  (1=1/2)  13C  nuclear  spins;  these  features  are  visible  at  correspondingly  lower 
contrast  in  both  (a)  and  (b). 


analytically  solving  7-Lys  (see  App.  B).  Double-electron  spin  flip  anti-crossings  are 
centered  at  Ep  =  0,  ±A\\h/ ge/iB  (Fig.  4-2b)  with  anti-crossing  coupling  strengths 
corresponding  to  those  predicted  by  simulation  and  plotted  in  Fig.  4-2a.  Note  the 
excellent  agreement  between  experiment  and  theoretical  model. 


63 


Figure  4-5:  Detail  of  double-electron  and  electron-nuclear  spin  flip  anti-crossings  at 
B  l  =  30  G,  where  greyscale  represents  normalized  fluorescence.  The  RF  excitation 
power  was  reduced  by  10  dB  to  mitigate  power  broadening.  This  allows  the  higher- 
resolution  features  to  be  observed,  however,  it  is  at  the  expense  of  signal-to-noise 
ratio.  Electron-nuclear  spin  flip  anti-crossings  of  coupling  strength  Ng  ~  250  kHz  are 
seen  centered  around  Np  ~  ±0.35  G,  separated  in  frequency  by  Pgs  at  2.877  GHz  and 
2.872  GHz.  Numerical  solutions  to  B.gs  are  plotted  over  the  data,  with  the  diameter 
of  the  circles  correspond  to  the  probability  of  the  transition. 

4.5  Conclusions  and  Outlook 

4.5.1  Utility  for  sensing 

The  numerical  simulation  of  %gs  in  the  <  40  G  transverse-field  regime  gives  insight 
into  the  level  structure  and  the  region  where  the  crossings  and  anti-crossings  can  be 
seen;  here  we  describe  several  useful  applications  in  sensing  and  atomic  memories. 
First,  the  double-electron  spin  flip  anti-crossing  addresses  the  ever-present  challenge 
of  isolating  sensitivity  to  multiple  fields  by  suppressing  sensitivity  to  changes  in  axial 
magnetic  fields.  As  the  applied  transverse  magnetic  held  is  increased,  the  strength 
of  this  anti-crossing  increases  quadratically,  and  the  range  of  B\\  to  which  the  NV 
ensemble  is  insensitive  increases  [62],  This  is  especially  useful  for  electric  held  and 
temperature  sensing,  where  the  effects  of  magnetic  helds  can  limit  sensitivity. 

Secondly,  the  double-electron  spin  hip  anti-crossing  causes  two  of  the  three  hyper- 
hne  transition  levels  to  cross.  At  these  degeneracies,  ESR  spectra  experience  doubled 
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signal  contrast.  In  CW  diamond-based  sensing  applications,  where  sensitivity  is  lin¬ 
early  proportional  to  the  ESR  resonance  contrast,  twice  the  sensitivity  is  expected 
in  this  regime.  Furthermore,  as  this  anti-crossing  strength  increases  and  the  hyper- 
fine  transition  levels  are  compressed  in  energy  spacing,  all  three  nuclear  levels  can 
be  effectively  degenerate,  providing  up  to  a  factor  of  three  increase  in  sensitivity  to 
electric  fields  or  temperature.  Up  until  now,  we  have  shown  experimental-theory 
agreement  for  the  six  hyperhne  transitions  present  in  a  14N  sample,  where  the  nuclear 
spin  InN  =  1.  Note  that  the  four  transitions  present  in  15N  allow  for  a  complete 
nuclear  spin  degeneracy  at  zero  applied  axial  magnetic  field  (see  App.  A). 

Next,  we  propose  a  transverse-field  method  to  remove  temperature  sensitivity 
from  diamond-based  magnetometry  measurements.  Taking  the  derivative  of  uc  with 
respect  to  temperature,  the  temperature-dependent  Dgs  term  in  the  denominator  of 
the  second  order  term  leads  to  a  decrease  in  duc/dT  with  increased  transverse  mag¬ 
netic  held  [35].  We  expect  a  1%  decrease  in  duiJdT  with  80  G  applied  transverse 
magnetic  held.  Combining  this  advantage  with  the  aforementioned  increase  in  con¬ 
trast  and  insensitivity  to  changes  in  Bz  in  a  high  signal-to-noise  ratio  light-trapping 
diamond  waveguide  results  in  a  regime  that  is  well-suited  for  temperature-stabilized 
measurements,  relying  on  the  ability  to  probe  multiple  sub-ensembles  with  different 
crystallographic  projections  of  the  applied  transverse  magnetic  held.  Using  these 
multiple  sub-ensembles  with  known  temperature  dependences  allows  the  tempera¬ 
ture  dependence  to  be  stabilized,  similar  to  the  method  proposed  by  Hodges  et  al. 
using  strain  engineering  [22],  This  method  further  decouples  the  diamond  from  its 
thermal  environment,  addressing  an  ongoing  challenge  in  diamond-based  electric  and 
magnetic  held  sensing. 

4.5.2  Implementation  with  15N 

Unlike  14N,  which  is  a  spin-1  system,  15N  is  a  spin-1/2  system.  The  hyperhne  splitting 
of  the  NV  center  would  therefore  result  in  two  resonances  (±1/2)  instead  of  three  res¬ 
onances  (0,  ±1).  This  gives  a  50%  increase  in  contrast,  which  is  linearly  proportional 
to  sensitivity.  Furthermore,  at  zero  applied  axial  magnetic  held  (B ||  =  0)  full  contrast 
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is  achieved,  even  at  non-zero  applied  transverse  magnetic  fields.  The  electron-nuclear 
spin  flip  anti-crossing  is  not  present  in  the  15N  case.  Under  a  transverse  magnetic 
field,  the  two  nuclear  spin  states  mix.  Therefore  the  application  of  a  linearly  polarized 
microwave  field  induces  transitions  between  all  (2  x  2)  possible  states,  giving  rise  to 
the  doublets  of  the  ODMR  spectrum  in  Fig.  4-6a. 
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Figure  4-6:  Comparison  of  15N  (a)  and  14N  (b)  samples  at  B±  =  30  G. 


4.5.3  Adiabatic  passage  protocol  for  nuclear  spin  polarization 

Using  these  low-transverse-field  anti- crossings,  we  introduce  an  alternate  route  to 
polarize  the  nuclear  spin  host  of  the  nitrogen  vacancy  center  [72,  75,  76,  77,  78]. 
This  method  articulates  successive  shifting  of  the  magnetic  field  and  radio-frequency 
pulses  iteratively  to  transfer  the  spin  polarization  into  a  single  nuclear  state,  similar 
to  methods  explored  for  optical  quantum  dots  [79,  80,  81].  Unlike  existing  schemes 
at  ~  500  G  and  ~  1000  G  axial  magnetic  fields  for  the  excited  state  and  ground  state 
LACs,  this  approach  relies  on  control  of  weak  magnetic  fields  and  state-selective  RF 
pulses.  It  requires  fewer  pulses  and  a  smaller  range  of  RF  excitation  than  other  low- 
magnetic  field  recursively  repeated  protocols  [82],  The  heart  of  the  scheme  relies  on 
adiabatic  passage,  which  results  in  a  change  in  nuclear  spin  state.  Following  Landau- 
Zener-Stueckelberg  theory  [83,  84,  85],  the  probability  P  of  adiabatic  passage  through 
an  anti-crossing  with  coupling  strength  A  is  P  =  1  —  e^~nhA2  where  ^  is  the 
rate  of  change  of  energy  difference  as  the  gap  is  approached.  Under  these  conditions, 
sweeping  a  bias  field  through  the  ±B\\  =  Np  electron- nuclear  spin  flip  anti-crossing 
results  in  spin  exchange  from  \ms  =  =pl ,  m/  =  ±1)  to  \ms  =  ±l,mj  =  0).  For  a 
99.99%  probability  of  adiabatic  passage,  a  rate  of  magnetic  field  change  of  1  G/225 
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/xs  is  required. 

While  many  nuclear  spin  exchange  sequences  are  possible,  here  we  note  a  two- 
step  example  to  polarize  into  mi  =  0.  Assuming  an  initial  ensemble  of  NVs  in  a 
superposition  of  eigenstates  1 0,  —1),  1 1, 1),  and  1 1, 0) ,  which  can  be  prepared  through 
state-selective  excitations,  a  single  application  of  the  protocol  results  in  a  superpo¬ 
sition  of  eigenstates  1 0,  —  1) ,  1 1, 0) ,  and  |  — 1,0).  The  population  transferred  through 
the  electron-nuclear  spin  flip  anti-crossing  is  then  shelved  in  ms  =  0  and  the  remain¬ 
ing  population  in  |0,  — 1)  is  transferred  to  |  —  1 ,  —  1)  with  a  resonant  RF  pulse.  The 
second  application  of  the  protocol  then  results  in  eigenstates  | — 1, 0),  1 0,  0),  and  1 1,  0) . 
Exploiting  the  robustness  of  the  nuclear  spin  state  against  optical  excitation  of  the 
NV,  all  ms  levels  can  be  pumped  into  the  ms  =  0  with  green  laser  excitation  [86]. 
This  type  of  protocol  would  be  useful  both  for  environmental  sensing  applications, 
due  to  the  increased  contrast,  but  also  for  applications  like  diamond-based  gyroscopes 
[40,  87]  and  atomic  memories  [88,  78,  89],  which  utilize  the  longer-lived  nuclear  spin 
state. 

4.5.4  Next  steps  and  future  applications 

The  transverse  magnetic  fields  may  help  overcome  the  influence  of  inhomogeneous 
strain  present  from  diamond  growth  processes,  a  challenge  in  working  with  very  large 
ensembles  of  NV  centers.  Combining  this  with  insensitivity  to  small  axial  magnetic 
fields  and  potential  insensitivity  to  temperature,  the  resultant  system  also  presents 
itself  as  a  strong  candidate  for  a  quantum  memory.  Transverse  and  axial  magnetic 
fields  could  be  used  to  read  and  write  to  the  system,  which  is  otherwise  isolated  from 
its  environment. 

In  conclusion,  in  this  work  we  present  a  study  of  the  effect  of  weak  transverse 
magnetic  fields  on  the  NV  system.  We  predict  and  experimentally  verify  an  electron- 
nuclear  spin  flip  anti-crossing,  with  close  agreement  between  theory  and  experiment. 
We  also  experimentally  measure  the  predicted  double-electron  spin  flip  anti-crossing 
under  transverse  magnetic  fields.  These  anti-crossings  show  potential  for  a  variety 
of  sensing  applications  because  of  (i)  increased  signal  contrast,  (ii)  insensitivity  to 
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axial  magnetic  fields  near  the  anti-crossing  points,  (iii)  potential  for  nuclear  spin 
polarization  schemes,  and  (iv)  selective  decoupling  of  the  NV  from  its  environment 
through  the  use  of  different  orientations  with  varying  transverse  fields.  With  applied 
transverse  magnetic  fields  on  the  order  of  tens  of  gauss,  we  can  achieve  increases  in 
contrast  on  the  order  of  those  seen  previously  only  at  the  much  higher,  axially-applied 
~  500  G  and  ~  1000  G  LACs. 
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Chapter  5 


Towards  a  diamond-based  frequency 
standard 

5.1  Introduction 

Precision  clocks  and  oscillators  are  crucial  for  a  range  of  applications  including  nav¬ 
igation,  computing  systems,  and  communications.  Some  of  the  most  accurate  time 
and  frequency  standards  developed  to  date  are  atomic  clocks,  which  derive  their  sta¬ 
bility  from  electronic  transitions  in  vapor-phase  atoms.  Most  modern  atomic  clocks 
rely  on  precisely  measured  internal  hyperhne  level  splitting  of  atoms.  However,  as  the 
focus  is  optimizing  the  stability,  such  clocks  often  require  a  volumes  of  up  to  several 
tens  of  cubic  meters,  too  large  for  many  applications.  At  the  other  extreme,  research 
has  been  directed  toward  miniaturizing  atomic  clocks  for  applications  such  as  satel¬ 
lites  and  in  situ  scientific  instruments.  Much  research  has  focused  on  miniaturizing 
such  systems,  resulting  in  chip-scale  vapor-cells  [90],  ultra-cold  atomic  traps  [91],  and 
vertical-cavity  surface-emitting  lasers  [92],  Moreover,  commercially  available  compact 
cesium  clocks  now  achieve  a  fractional  frequency  deviation  of  6  x  10_10/\/Hz[93], 
representing  roughly  100  times  better  timing  stability  than  temperature-stabilized 
quartz-based  oscillators  in  a  package  only  17  cc  in  volume. 

Concerted  research  efforts  to  miniaturize  atomic  clocks  have  resulted  in  on-chip  de¬ 
vices  employing  small  vapor  cells  that  are  integrated  with  detectors  and  laser  sources 
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[94,  95].  Solid  state  alternatives  could  not  only  reduce  the  volume,  but  could  also 
enable  better  performance  in  harsh  environments.  Smaller  devices  could  enable  high- 
performance  clock  generators  for  short-range  next-generation  wireless  (millimeter- 
wave)  communications,  since  existing  quartz-based  oscillators  suffer  from  high  phase 
noise  at  high  frequencies  [96].  In  addition,  mobile  devices  can  derive  an  accurate 
time  signal  from  GPS  for  maintaining  communication  standards  [97],  but  are  sus¬ 
ceptible  to  obstruction  or  jamming  of  the  external  clock  signal  —  again,  a  portable 
high-performance  clock  that  is  small  enough  to  be  incorporated  into  mobile  devices 
could  maintain  synchronization. 


A  recent  proposal  suggests  a  pulsed  protocol  that,  addressing  a  very  large  en¬ 
semble  of  NV  centers,  could  result  in  a  solid-state  clock  approaching  the  stability  of 
chip-scale  atomic  clocks  [22],  This  protocol  for  a  stable  oscillator  based  on  magnetic 
sub-level  transitions  in  the  diamond  NV  center.  This  protocol  relies  on  exploiting  the 
NV’s  long  spin  coherence  time  under  rephasing  sequences  and  using  a  large  number  of 
NV  centers.  Due  to  the  NV’s  long  spin  coherence  time,  a  high  density  of  spins  in  the 
solid,  and  optical  spin  detection,  it  is  estimated  that  its  frequency  stability  could  rival 
the  performance  of  the  newest  chip-scale  Cs  and  Rb  standards  in  a  package  that  is 
at  least  2  orders  of  magnitude  smaller  and  lighter.  Developing  an  atom-like  standard 
in  a  solid-state  host  could  enable  its  integration  into  semi-conductor  fabrication  pro¬ 
cesses,  which  could  lead  to  rapid  reductions  in  volume  and  power  requirements.  The 
light-trapping  diamond  waveguide  (LTDW)  is  an  excellent  candidate  for  addressing 
large  ensembles  of  NV  centers.  The  proposed  clock  protocol  is  a  spin-echo  technique 
that  eliminates  magnetic  held  effects  to  first  order,  but  is  sensitive  to  fluctuations  in 
temperature.  In  this  chapter,  several  methods  to  combat  these  temperature  fluctua¬ 
tions  are  described,  including  strain  engineering  and  the  application  of  a  transverse 
magnetic  held.  The  feasibility  of  practical  diamond-based  timekeeping  is  assessed, 
and  the  steps  towards  its  implementation  are  outlined. 
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5.2  Timekeeping  protocol 


The  protocol  relies  on  the  manipulation  of  the  spin-1  NV  ground  state  [98],  as  il¬ 
lustrated  in  Figure  1-1.  The  NV  center  is  first  prepared  into  the  spin-1  ground 
state  |0).  A  pulsed  technique  is  used  to  monitor  the  frequency  drift.  Similar  to 
the  Ramsey  spectroscopy  used  in  atomic  clocks,  this  results  in  a  fluorescence  re¬ 
sponse  that  varies  with  respect  to  frequency  drift,  8uj.  The  clock  protocol  mod¬ 
ifies  a  traditional  spin-echo  protocol,  changing  the  flip-angle  of  the  first  and  last 
pulses  from  7t/2  to  7t/4.  After  performing  the  clock  protocol,  the  NV  is  in  state 
I'lpf)  =  Uciock  IV’o)  =  sinfAarr)  |0)  —  |  cos(5o;r)(|+l)  +  |— 1)).  This  drift-dependent 
sinusoidal  signal  acts  as  a  long-time-scale  correction  to  a  lower  stability  voltage  con¬ 
trolled  oscillator  (VCO)  at  ~  3  GHz. 

The  figure  of  merit  for  the  stability  of  a  frequency  standard  is  fractional  frequency 
deviation.  This  is  a  measure  of  how  much  an  oscillator  drifts  around  the  oscillation 
frequency.  Equation  5.1  describes  how  the  fractional  frequency  deviation  is  related 
to  the  observed  fluorescence  signal  (M). 


5cu  _  1  AM 
w0  w0  |  d{M)/du 


(5.1) 


where  AM  =  J (M2)  —  (M)2  is  the  standard  deviation  of  the  operator  describing 
the  spin  expectation  value  for  the  fluorescence  measurement  and  \d(M)/du\  is  the 
maximum  slope  of  the  the  signal  as  it  changes  with  frequency  drift. 

To  facilitate  comparison  between  systems  and  recognize  that  multiple  measure¬ 
ments  will  improve  the  SNR  term,  we  can  simplify  this  to  Equation  5.2,  where  r  is 
time  separation  between  the  pulses,  T  =  Mr  is  the  total  interrogation  time  for  M 
measurements,  and  N  is  the  number  of  NVs  in  the  ensemble.  Similar  to  magnetic  field 
sensitivity,  maximizing  the  collection  efficiency  of  the  system  as  well  as  the  contrast 
will  improve  this  metric. 


5u>\ 

w0  / m,n  DqsVtNT 


(5.2) 
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5.3  Demonstration 


5.3.1  Single  NV  clock 


delay  (ns) 


Figure  5-1:  a)Raster  scan  of  the  fluorescence  of  the  diamond  sample  when  pumped 
with  532  nm  green  light.  The  NV  used  is  noted  in  red.  b)  g ^  plot  for  this  NV, 
verifying  that  it  is  a  single  defect,  c)  CW  ESR  plot  for  this  defect  center,  indicating 
that  the  resonance  frequency  is  at  2.87  GHz  as  is  typical,  d)  Diagram  of  setup  used 
to  perform  measurements  on  this  defect  center. 

Here,  we  demonstrate  the  NV  protocol  on  a  single  NV  center,  demonstrating 
an  fractional  frequency  deviation  of  2  x  10~6/\/Hz-  To  confirm  that  the  NV  being 
interrogated  is  a  single  emitter,  a  fluorescence  image  (Fig.  5-la),  a  scan  (Fig. 
5-lb),  and  aCW  ESR  scan  (Fig.  5-lc)  are  recorded.  The  experimental  setup  is  shown 
in  Fig.  5-ld.  Using  a  confocal  microscope,  we  excite  the  NV  center  using  a  532  nm 
green  laser.  Microwave  pulses  are  applied  using  a  15  /mi  wire  approximately  20  /mi 
from  the  target  NV.  The  photoluminescence  signal  is  spectrally  filtered  from  the  input 
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signal  and  sent  to  the  detector.  A  slight  splitting  in  the  CW  ESR  signal  for  this  NV 
center  is  observed,  which  may  be  caused  by  a  small  magnetic  held.  To  first  order, 
however,  this  held  should  not  affect  the  operation  of  the  clock  protocol.  This  slight 
magnetic  held  shows  up  as  the  beating  frequencies  visible  in  the  plot  of  the  Rabi 
oscillations  (Figure  5-2).  From  this  plot,  we  chose  the  7t/4  time  to  be  22  ns  and  the 
7 r  time  to  be  88  ns. 


Figure  5-2:  Rabi  oscillations  of  single  NV  center 


Substituting  directly  into  Equation  5.1,  using  the  standard  deviation  of  the  7t/4 
point  on  the  curve  and  the  slope  around  it,  we  find  a  fractional  frequency  deviation 
of  2  x  10  6  for  r  =  10  /is.  This  corresponds  to  a  thermal  stability  of  ~  80  mK,  which 
is  reasonable  for  a  diamond  sample  with  no  thermal  management,  and  indicates  that 
the  single  NV  case  is  likely  limited  by  thermal  huctuations. 


5.3.2  Small  ensemble  clock 

The  hrst  step  towards  the  demonstration  of  the  clock  protocol  on  an  NV  ensemble 
is  the  implementation  of  pulsed  ESR  and  Rabi  oscillations.  Here,  a  small  ensemble 
(~  106  NVs)  is  implemented  in  a  non-LTDW  geometry.  Using  a  single-pass,  confocal 
setup,  clear  Rabi  oscillations  are  seen. 
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Figure  5-3:  The  clock  protocol  is  performed  on  a  single  NV  center  in  diamond.  The 
fits  have  r-square  values  of  0.73  and  0.8  for  r  =  8  /xs  (blue)  and  r  =  10  /xs  (red) 
respectively. 

In  Figure  5-4,  Rabi  oscillations  are  shown  on  the  left  and  the  clock  protocol  is 
shown  on  the  right,  for  r  =  2  /xs  (blue)  and  r  =  4  /xs  (red).  For  the  fits,  the  height, 
contrast,  and  location  of  zero  detuning  were  adjusted,  but  the  frequency  of  oscillation 
was  not.  The  fractional  frequency  deviation  is  calculated  to  be  ~  10-7,  a  minor 
improvement  over  single  NV  clock.  This  is  consistent  with  both  implementations 
being  limited  by  ambient  temperature  fluctuations. 


5.4  Pulsed  protocols  using  a  LTDW 

5.4.1  RF  Uniformity  &  Antenna  development 

Non-uniformity  of  the  applied  RF  field  in  pulsed  protocols  results  in  a  degradation 
of  the  signal,  since  RF  power  affects  the  pi  pulse  time.  This  is  particularly  relevant 
in  systems  like  the  LTDW  where  the  NVs  are  distributed  over  a  large  area.  For 
ensembles  on  the  order  of  100  /xm2,  a  mm-scale  loop  of  wire  provides  adequate  im- 
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(a)  (b) 


Figure  5-4:  (a)  Demonstration  of  Rabi  oscillations,  (b)  Demonstration  of  the  clock 
protocol  on  small  ensemble  of  NV  centers,  with  the  y-axis  shows  normalized  fluores¬ 
cence  and  x-axis  shows  detuning  from  resonance  in  Hz 


provement  over  the  straight  length  of  wire  commonly  used  in  single-NV  experiments. 
This,  however,  is  inadequate  for  NVs  spread  out  across  several  millimeters,  as  in  the 
case  of  the  LTDW.  A  more  uniform  held  can  be  realized  with  a  Helmholtz  coil  geom¬ 
etry.  However,  the  inductance  of  the  coils  prevents  GHz  held  switching.  Recent  work 
has  shown  spatial  uniformity  across  millimeter-scale  areas  [99,  100]  using  split  ring 
resonator  antennas.  These  antennas  are  extremely  narrow  band  and  not  suited  for 
magnetometry  experiments  requiring  an  applied  bias  held.  Other  designs  use  more 
broadband  "lollipop"  type  geometries  to  accommodate  high-dynamic-range  measure¬ 
ments  [101].  Adequate  power  delivery  remains  a  challenge  for  large-ensemble-based 
NV  sensing.  Additionally,  alternative  methods  of  exciting  the  transitions  have  been 
explored  [102],  Lab-based  investigation  into  alternative  antenna  designs,  ranging  from 
PCB-based  loops  and  spirals  to  coax-loops  with  a  variety  of  grounding  configurations 
to  three  dimensional  resonant  cavities  have  not,  to  date,  yielded  a  result  with  signif¬ 
icant  improvement  over  a  4  mm  loop  of  wire.  As  such,  the  demonstrations  of  pulsed 
protocols  in  Section  5.4.2  are  extremely  limited  by  the  non-uniform  RF  power  across 
the  sample. 

One  method  to  overcome  the  lack  of  uniformity  without  engineering  a  custom 
antenna  is  to  use  advanced  dynamical  decoupling  [103]  and  optimal  control  techniques 
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[104],  such  as  those  designed  using  the  DYNAMO  package  [105,  106].  We  estimate 
that  these  techniques  are  most  helpful  if  the  RF  inhomogeneity  is  limited  to  a  few 
percent  across  the  NV  volume,  otherwise  the  Rabi  oscillation  signal  from  one  part  of 
the  sample  can  directly  cancel  that  from  another  part  which  experiences  a  180  degree 
phase  shifted  pulse. 

5.4.2  Demonstration  of  Rabi  oscillations,  Hahn  Echo 

Given  the  known  limitations  due  to  the  RF-inhomogeneity  of  a  4  mm  loop  antenna 
used  to  apply  RF  to  the  3  mm  square  LTDW  sample,  it  is  not  surprising  that  the  re¬ 
sulting  Rabi  oscillations  from  the  ensemble  begin  to  cancel  each  other  rapidly.  Figure 
5-5  shows  Rabi  oscillations  as  the  microwave  pulse  is  increased  from  10  ns  to  3000  ns. 
This  is  the  first  step  towards  determining  the  correct  duration  for  7 r  pulses  and  other 
pulses  as  the  basis  for  more  complex  sensing  applications. 
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Figure  5-5:  Rabi  sequence  performed  on  the  LTDW  platform 

Next,  a  pulsed  ODMR  experiment  is  implemented.  Taking  the  7 r  rotation  time 
from  the  Rabi  oscillation  plot  to  achieve  maximum  contrast  (~  710  ns)  we  sweep  the 
microwave  excitation  at  zero  applied  held.  Imperfect  RF  uniformity  results  in  visible 
sidebands. 

Figure  5-6  shows  a  Hahn  Echo  sequence  (7t/2  -  t  -  tt  -  t  -  7t/2)  with  revivals 
out  to  r  ~  200  microseconds  using  the  LTDW  sample.  The  spacing  of  the  revivals 
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Figure  5-6:  Hahn  echo  sequence  performed  on  LTDW  platform 


corresponds  to  the  total  applied  magnetic  field  of  roughly  35  G.  The  coherence  time 
measured  here  is  of  the  same  order  of  magnitude  as  coherence  times  measured  on  single 
nitrogen-vacancy  defect  centers  in  a  comparable  sample  on  a  confocal  microscope 
setup  previously.  However,  the  SNR  seen  in  this  result  is  low,  indicating  that  likely 
only  a  small  percentage  of  NVs  are  being  addressed  in  this  scan. 


5.5  Decoupling  from  temperature  effects 

5.5.1  Theoretical  improvement 

By  engineering  two  ensembles  of  NV  centers  where  each  zero  field  splitting  frequency 
has  a  different  different  dependence  on  temperature  =  /3n),  we  can  obtain 

a  fractional  frequency  deviation  that  is  a  factor  of  £  =  ^Jl  +  2(^/3l/ji  )2  worse  than 
that  of  the  ideal  (zero  temperature  fluctuation)  case.  Even  in  the  single  NV  clock, 
this  stability  limit  is  an  improvement  over  the  temperature-limited  stability.  To  get 
as  close  as  possible  to  the  ideal  result,  the  difference  between  /3i  and  should  be 
maximized. 

The  tetrahedral  lattice  of  the  diamond  provides  four  convenient  sub-ensembles 
of  NVs  that  are  experiencing  the  exact  same  temperature  fluctuations,  due  to  the 
excellent  thermal  conductivity  of  diamond.  Two  methods  of  altering  dDgs/dT  will  be 
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explored  in  this  section. 


5.5.2  Change  in  dDgs/dT  using  transverse  magnetic  fields 

By  solving  for  the  energy  eigenvalues  of  the  Hamiltonian  detailed  in  Section  1.1.2 
and  applying  second  order  energy  corrections,  we  find  that  the  zero  held  splitting 
frequency  can  be  described  by  ujzfs  =  Dgs  +  x  • 


Figure  5-7:  Change  in  temperature  dependence  with  the  application  of  transverse 
magnetic  held. 

The  temperature  dependence  of  one  sub-ensemble  of  NV  centers’  zero  held  split¬ 
ting  frequency,  /3  ~  —74  kHz/K,  can  be  altered  with  the  application  of  a  strong 
magnetic  held  perpendicular  to  the  NV  ensemble  axis  as  described  in  Eqn.  5.3  and 
plotted  in  Fig.  5-7. 


Pb± 


dDgs  f  ,  _  3/i0g  2  \  _  dPgs 
dT  V  2 h2Dgs  ±)  ~  dT 


(5.3) 


This  temperature  dependence  falls  oh  quadratically  with  increased  applied  trans¬ 
verse  magnetic  held.  A  1%  difference  between  /3B±. i  and  (3b±, 2  is  expected  with  the 
application  of  an  80  G  transverse  magnetic  held,  while  a  2%  difference  is  expected 
with  the  application  of  a  120  G  transverse  magnetic  held. 

B±  can  be  applied  uniformly  with  a  combination  of  Helmholtz  coils  and  a  Halbach 
array  of  permanent  magnets.  Resolving  a  10  Kelvin  temperature  shift  will  require 
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resolution  of  the  zero-field  splitting  frequency  to  ~  750  Hz  to  adequately  differentiate 
between  the  sub-ensemble  that  is  most  affected  by  B±  and  the  ones  that  are  minimally 
affected.  Once  a  sample  has  been  characterized,  measuring  the  difference  between 
the  zero  held  splitting  frequency  dependence  over  a  temperature  range  will  yield  an 
absolute  temperature  measurement,  immune  to  drift.  However,  since  the  strain  profile 
of  different  diamond  samples  vary,  each  device  will  need  to  be  initially  calibrated. 

5.5.3  Change  in  dDgs/dT  using  strain 

Anisotropically-applied  strain  can  also  be  used  to  change  the  temperature  depen¬ 
dence  of  one  or  more  of  the  diamond  NV  sub-ensembles  as  detailed  in  [22] .  While 
the  lack  of  a  large  permanent  magnetic  held  makes  this  is  an  attractive  option,  par¬ 
ticularly  for  a  compact  system  for  in  situ  measurements,  extremely  high  precision 
machining  is  required  to  produce  a  clamp  that  could  apply  sufficiently  uniform  strain 
across  the  sample.  Alternately,  carefully  engineered  diamond  growth  can  produce 
highly  strained  regions  of  the  sample  that  could  be  individually  probed.  Strain  inho¬ 
mogeneity  would  have  the  effect  of  broadening  the  resonance  linewidths,  decreasing 
the  stability  of  the  clock,  so  large  areas  of  uniform  strain  are  needed.  Furthermore, 
due  to  diamond’s  high  Young’s  modulus  [107],  large  forces  are  required  to  produce 
this  change,  which  can  break  or  shatter  the  brittle  diamond  samples  if  not  applied 
evenly.  Previous  work  in  using  diamond  for  pressure  sensing  indicates  sensitivity  on 
the  MPa/\/Hz  scale  [67,  108]. 

5.6  Conclusions  and  outlook 

5.6.1  Feasibility  of  a  competitive  diamond-based  frequency 
standard 

To  date,  temperature  instability  has  been  the  limiting  factor  for  a  practical  diamond- 
based  frequency  standard  competitive  with  chip-scale  atomic  clocks.  Table  5.1  sum¬ 
marizes  temperature  stability  corresponding  to  fractional  frequency  deviation,  given 
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an  phenomenologically  derived  value  of  dDgs/dT  of  75  kHz/K  and  an  average  value 
of  Dgs  ~  2.87  GHz. 


Table  5.1:  Required  temperature  stability 


Temp.  Stability 

Fractional  Frequency  Deviation 

1  K 

~  2.6  x  10"5 

1  mK 

~  2.6  x  10"8 

1  /jK 

~  2.6  x  IQ’11 

For  GPS,  a  fractional  frequency  deviation  of  better  than  ~  10-12  is  needed,  and 
better  than  ~  1CP14  is  needed  for  spacecraft  navigation.  For  comparison,  the  frac¬ 
tional  frequency  deviations  of  several  types  of  atom-based  clocks  are  given  in  Table 
5.2. 


Table  5.2:  Comparison  with  atom-based  frequency  standards 


Atom  type 

Fractional  Frequency  Deviation 

Cs 

~  io-13 

Rb 

~  10“15 

H 

~  10“15 

St 

~  10-17 

Two  methods  have  been  explored  to  improve  the  stability  of  a  diamond-based 
frequency  standard  earlier  in  the  chapter.  First,  by  engineering  two  ensembles  of  NV 
centers  where  each  zero  held  splitting  frequency  has  a  different  different  dependence 
on  temperature  ,  we  can  obtain  a  fractional  frequency  deviation  that 

is  a  factor  of  £  worse  than  that  of  the  ideal  (perfect  temperature  stability)  case, 
where  f  =  ■  A  1%  difference  between  ft  and  ft  is  expected  with  the 

application  of  an  80  G  transverse  magnetic  held,  while  a  2%  difference  is  expected 
with  the  application  of  a  120  G  transverse  magnetic  held. 

This  ideal  fractional  frequency  deviation  for  an  ensemble  can  be  extrapolated  from 
the  non-temperature-limited  single  NV  result  of  ^  ~  10~6.  Improvement  scales  as 
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Table  5.3:  Factor  away  from  ideal  fractional  frequency  deviation 


P2-P1 

Pi 

£ 

0.01 

141 

0.02 

70 

0.03 

47 

0.10 

14 

y/N  where  N  is  the  number  of  NVs.  For  a  fractional  frequency  deviation  of  ~  1CT12, 
given  a  l%-2%  difference  in  temperature  dependence  of  the  zero  held  splitting  for 
two  same-temperature  ensembles,  an  ensemble  of  1016  ensembles  is  required.  This 
corresponds  to  a  density  of  33  ppm  in  a  3  mm  x  3  mm  x  0.3  mm  sample,  which 
is  not  commercially  available  at  this  time.  Currently  available  samples  at  densities 
approaching  this  have  low  quality  NVs  due  to  lattice  damage  and  other  impurities. 
Furthermore,  at  this  high  density,  the  higher  absorption  would  limit  the  path  length 
in  an  LTDW  geometry,  making  an  even  smaller,  denser  sample  preferable. 

Considering  the  challenges  in  implementation  of  pulsed  protocols  on  the  LTDW 
platform  due  to  RF  uniformity,  it’s  also  worth  considering  the  CW  clock  case.  By 
studying  the  Zeeman-split  ODMR  spectra  and  monitoring  the  frequency  of  each  of 
these  resonances,  two  signals  emerge,  as  shown  in  Section  3.2.2.  The  sum  of  the 
two  frequencies  gives  a  temperature  dependent  signal  that,  as  a  clock,  is  limited  as 
indicated  in  Table  5.1.  The  difference  of  the  two  frequencies  gives  a  temperature- 
independent  signal.  Given  previous  LTDW  magnetometry  results  [39],  temperature 
insensitivity  corresponding  to  a  fractional  frequency  deviation  of  ~  10-9  has  been 
achieved.  However  this  is  not  suitable  for  a  frequency  standard  unless  it  can  be 
magnetically  shielded  extremely  well,  as  it  is  sensitive  to  nT  level  magnetic  fields, 
unlike  the  clock  protocol,  which  is  insensitive  to  magnetic  fields  to  first  order. 

In  conclusion,  there  are  three  hurdles  that  must  be  overcome  in  order  to  implement 
a  diamond-based  frequency  standard  that  is  competitive  with  compact  atomic  clocks: 
i)  uniformly  applied  RF,  ii)  decoupling  from  ambient  temperature  fluctuations,  and 
iii)  production  of  an  adequately  high  NV  density  sample  with  good  quality  NVs. 
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Several  methods  for  combatting  the  first  two  issues  have  been  discussed,  and  with 
the  rapid  advances  in  diamond-growth  processes  seen  over  the  past  several  years,  its 
not  unthinkable  to  imagine  that  a  competitive  diamond-based  frequency  standard 
could  be  seen  in  the  next  several  years. 
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Chapter  6 


Conclusions  and  outlook 


After  a  brief  introduction,  this  thesis  motivated  and  presented  the  light-trapping  di¬ 
amond  waveguide  device,  and  demonstrated  its  utility  for  addressing  large  ensembles 
of  NV  centers.  This  simple  device’s  exceptional  excitation  efficiency  and  good  signal 
collection  enables  in  excess  of  5%  conversion  efficiency  of  pump  photons  into  opti¬ 
cally  detected  magnetic  resonance  (ODMR)  fluorescence,  a  three  orders  of  magnitude 
improvement  over  previous  single-pass  geometries.  This  dramatic  enhancement  of 
ODMR  signal  enables  precision  broadband  separable  measurements  of  magnetic  field 
and  temperature  in  the  low-frequency  range  inaccessible  by  dynamical  decoupling 
techniques.  The  increased  sensitivity  achievable  through  ensemble-based  magnetic 
held  sensing  lends  itself  to  many  existing,  mature  applications,  where  single-NV  mag¬ 
netic  held  sensing  remains  limited  to  the  lab  setting.  In  particular,  geoscience  and  ob¬ 
ject  detection  are  two  application  spaces  whereas  ensemble-based  sensing  shows  great 
potential.  Others  include  biomagnetism,  bioassays,  condensed  matter  physics,  mate¬ 
rial  science,  biothermometry,  bulk  magnetometry  for  surveying,  and  hyper-polarized 
media  for  NMR 

The  LTDW  is  then  used  to  study  of  the  effect  of  weak  transverse  magnetic  helds 
on  the  NV  system.  An  electron- nuclear  spin  hip  anti-crossing  is  predicted  and  ex¬ 
perimentally  verihed.  The  previously  predicted  double-electron  spin  hip  anti-crossing 
under  transverse  magnetic  helds  is  also  measured.  These  anti-crossings  show  poten¬ 
tial  for  a  variety  of  sensing  applications  because  of  (i)  increased  signal  contrast,  (ii) 
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insensitivity  to  axial  magnetic  fields  near  the  anti-crossing  points,  (iii)  potential  for 
nuclear  spin  polarization  schemes,  and  (iv)  selective  decoupling  of  the  NV  from  its 
environment  through  the  use  of  different  orientations  with  varying  transverse  fields. 

This  sets  the  stage  for  a  large-ensemble  device  that  could  be  used  for  a  realistic 
diamond-based  clock,  pending  advances  in  diamond  material  growth  and  processing. 
The  proposed  clock  protocol  is  demonstrated  on  a  single  NV  and  on  a  small  ensemble 
of  NVs,  followed  by  a  discussion  of  the  remaining  challenges  towards  implementation 
using  the  LTDW  or  other  large-ensemble  platform.  The  two  main  challenges  can 
be  summarized  as  difficulties  in  obtaining  applied  held  uniformity  and  decoupling 
from  the  environment.  Beyond  the  simple  method  of  removing  temperature  effects 
by  monitoring  both  halves  of  the  split  ODMR  signal,  two  more  advanced  methods 
for  combatting  ambient  temperature  fluctuations  are  detailed.  Finally,  the  feasibility 
of  a  competitive  diamond-based  frequency  standard  is  considered. 

One  limit  to  the  current  light-trapping  diamond  waveguide  setup  is  non-optimal 
fluorescence  collection  (around  20%).  By  collecting  light  from  more  sides  of  the 
diamond,  and  closer  to  the  diamond,  this  number  can  be  improved.  An  improvement 
in  sensitivity  of  \/N  for  an  V-fold  increase  in  the  number  of  photons  collected  is 
expected  due  to  a  decrease  in  the  photon  shot  noise  limit.  Future  iterations  of  the 
LTDW  device  might  incorporate  fiber  optic  bundles  and  light  pipes  that  can  be  placed 
near  the  diamond  to  collect  the  signal.  These  will  also  serve  to  reduce  electrical  and 
magnetic  noise  by  further  removing  the  photodetector  from  the  setup.  In  addition, 
strategically  deposited  DBR  mirrors  could  also  increase  the  fluorescence  collection 
and  device  sensitivity. 

In  many  ways,  the  LTDW  device  geometry  is  a  work-around  to  enable  sensing 
using  large  NV  ensembles  while  the  material  development  processes  mature.  The  size 
of  the  device  makes  it  particularly  susceptible  to  inhomogeneities  in  the  applied  RF 
excitation,  laser  excitation,  and  magnetic  bias  held.  The  comparatively  large  volume 
limits  the  spatial  resolution  of  the  device  as  a  sensor.  An  ideal  large-NV-ensemble 
sensor  platform  would  be  based  on  a  high  density  sample  with  excellent  defect  quality. 
Diamond  growth  using  isotropically  pure  12C  as  well  improvements  in  growth  speed 
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and  crystal  quality  could  make  such  a  sample  a  reality.  For  the  diamond  clock,  a  33 
ppm  sample  is  needed  based  on  the  current  size  of  sample.  An  ideal  clock  sample 
would  be  smaller  with  higher  density  to  achieve  the  same  number  of  defects,  or,  lower 
density  with  improved  quality  NVs  with  longer  coherence  times. 

Diamond-based  sensing  has  reached  a  point  in  the  arc  of  it’s  popularity  where 
many  are  wondering  whether  it  will  live  up  to  the  hype.  In  order  for  the  multitude  of 
lab-based  experiments  to  gain  a  stronghold  in  the  market  and  compete  with  existing 
commercial  sensors,  diamond-based  sensors  need  to  hit  sensitivity  thresholds  that 
match  or  exceed  commercially  available  technologies,  rather  than  focusing  on  new 
records  compared  with  other  diamond-based  devices.  There  must  be  a  focus  on 
designing  robust,  reliable  devices  that  lend  themselves  to  mass  manufacturing,  rather 
than  one-off  experiments. 

NV  centers  in  diamond  are  sensitive  to  many  environmental  parameters.  This 
presents  one  of  the  strongest  arguments  against  pursuing  the  NV  over  searching  for 
other  promising  defect  centers,  as  the  ideal  magnetometer  is  sensitive  to  magnetic 
fields  and  nothing  else.  However,  this  multi-factor  sensitivity  forms  the  basis  for  the 
possibility  of  a  multi-purpose  sensor  that  could  be  extremely  compact  and  robust. 
Similar  to  the  idea  of  using  two  sub-ensembles  as  co-located  clocks  to  combat  noise 
from  thermal  fluctuations  presented  in  chapter  5,  each  sub-ensemble  could  be  used  in¬ 
dependently.  Alternatively,  three  orientations  could  be  used  for  vector  magnetometry 
while  a  fourth  records  temperature  data. 

This  environmental  sensitivity  is  a  major  concern  for  the  diamond  electron-spin 
clock.  A  realistic  implementation  that  is  competitive  with  existing  technologies  will 
rely  on  precise,  careful  engineering  and  concerted  effort  at  developing  a  high  quality, 
high  density  NV  sample.  If  adequate  density  is  not  achieved,  time  will  also  need  to 
be  spent  addressing  inhomogeneities  in  laser,  RF,  and  magnetic  bias  fields  that  result 
from  the  device’s  large  volume.  While  it  is  unlikely  that  a  hybrid  diamond  clock 
magnetometer  will  be  in  every  airplane  in  five  years,  the  pace  of  progress  in  diamond- 
based  sensing  over  the  past  five  years  indicates  exciting  advances  are  forthcoming. 
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Appendix  A 


Polymer  Photonics 


The  following,  tangentially  related  investigations  were  also  completed  over  the  course 
of  my  doctoral  work. 

Polymers  have  appealing  optical,  biochemical,  and  mechanical  qualities,  including 
broadband  transparency,  ease  of  functionalization,  and  biocompatibility.  However, 
their  low  refractive  indices  have  precluded  wavelength-scale  optical  confinement  and 
nanophotonic  applications  in  polymers.  Here,  we  propose  a  suspended  polymer  pho¬ 
tonic  crystal  (SPPC)  architecture  that  enables  the  implementation  of  nanophotonic 
structures  typically  limited  to  high-index  materials.  Filters  and  integrated  networks, 
as  well  as  trace  gas,  force,  and  displacement  sensing  are  explored. 

A.l  Nanophotonic  Filters  and  Integrated  Networks 
in  Flexible  2D  Polymer  Photonic  Crystals 

This  work  appeared  in  the  following  publication:  [109] 

A. 2  High  sensitivity  gas  sensor  based  on  high-Q  sus¬ 
pended  polymer  photonic  crystal  nanocavity 

This  work  appeared  in  the  following  publication:  [110] 
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A. 2.1  Introduction 


We  present  high-sensitivity,  multi-use  optical  gas  sensors  based  on  a  one-dimensional 
photonic  crystal  cavity.  These  devices  are  implemented  in  versatile,  flexible  polymer 
materials  which  swell  when  in  contact  with  a  target  gas,  causing  a  measurable  cavity 
length  change.  This  change  causes  a  shift  in  the  cavity  resonance,  allowing  precision 
measurements  of  gas  concentration.  We  demonstrate  suspended  polymer  nanocavity 
sensors  and  the  recovery  of  sensors  after  the  removal  of  stimulant  gas  from  the  sys¬ 
tem.  With  a  measured  quality  factor  exceeding  104,  we  show  measurements  of  gas 
concentration  as  low  as  600  parts  per  million  (ppm)  and  an  experimental  sensitivity 
of  10  ppm;  furthermore,  we  predict  detection  levels  in  the  parts-per-billion  range  for 
a  variety  of  gases. 

Sensing  dilute  gases  is  of  interest  in  a  variety  of  fields  [111].  Much  research  has  been 
directed  at  developing  optical  sensors  for  this  task,  as  they  enable  a  high  signal-to- 
noise  ratio,  can  operate  in  ambient  conditions  over  a  wide  dynamic  range,  are  immune 
to  electromagnetic  interference,  and  can  be  compact  and  lightweight  [112,  113].  In 
particular,  optical  resonators  can  dramatically  enhance  the  light-matter  interaction 
time,  which  can  lead  to  improved  sensitivity  to  dilute  gas  concentrations  [114,  115]. 
Previous  work  has  employed  optical  fields  confined  in  inorganic  materials,  interacting 
evanescently  with  gas-sensitive  coatings,  such  as  polymers  [116,  117,  118,  119,  120], 
or  directly  with  the  gas  [121]  as  a  result  of  refractive  index  change.  In  this  work, 
we  use  optical  modes  confined  in  polymer  photonic  crystals  which  swell  in  response 
to  gas  exposure,  as  illustrated  in  Figure  1.  Since  the  mode  is  fully  in  the  polymer, 
even  minute  length  or  refractive  index  changes  of  the  cavity  resulting  from  interaction 
with  a  gas  cause  strong  and  measurable  changes  in  the  cavity  resonance  frequency 
(see  Fig.  A-l). 

A. 2. 2  Device  design 

The  device  is  based  on  a  one  dimensional  (ID)  photonic  crystal  fabricated  in  poly(methyl 
methacrylate)  (PMMA),  which  is  an  inexpensive,  flexible,  and  resilient  long-chain 
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Figure  A-l:  The  measurement  setup  and  a  conceptual  schematic  of  the  operation  of 
this  gas  sensor.  As  the  polymer  photonic  crystal  cavity  swells,  its  resonance  peak 
shifts.  Various  surface  treatments  will  enhance  sensitivity  to  particular  reagents. 


polymer  [109].  The  swelling  of  the  material  is  caused  by  the  interaction  of  the  PMMA 
ester  pendant  groups  via  hydrogen  bonding  with  organic  solvents  such  as  alcohols.  To 
interact  with  other  gas  molecules,  PMMA  can  be  functionalized  appropriately  [122], 
Conveniently,  PMMA  is  also  an  electron  beam  lithography  (EBL)  resist,  making  the 
device  very  straightforward  to  fabricate;  it  requires  no  mask  transfers  or  extra  wet 
or  dry  etching  steps  to  process,  unlike  semiconductor  materials  that  are  traditionally 
used  for  nanophotonics.  However,  due  to  the  low  refractive  index  of  PMMA  (mea¬ 
sured  to  be  ~  1-52),  transverse  confinement  of  light  by  total  internal  reflection  (TIR) 
is  difficult  to  achieve  when  clad  in  another  material.  The  suspended  one-dimensional 
photonic  crystal  resonators  described  here  overcome  this  problem  by  surrounding  the 
region  of  interest  with  air  on  all  sides.  This  also  allows  the  material  to  distend  without 
obstruction  when  the  analyte  is  present.  Periodic  rectangular  holes  in  the  structure 
form  a  distributed-Bragg  reflector  (DBR).  Light  is  confined  transversely  through  TIR 


and  longitudinally  via  a  photonic  bandgap.  Similar  suspended  architectures  have  been 
described  for  various  materials  including  high-index  semiconductors  [123,  124,  125] 
and  silicon  dioxide  [126].  Fig.  A-2a  shows  the  polymer  photonic  crystal  lattice,  which 
exhibits  a  bandgap  when  0.4  <  c o  <  0.45,  where  u  is  the  normalized  frequency,  as  ob¬ 
tained  from  finite  difference  time  domain  (FDTD)  simulations  [109].  A  defect  region 
creates  a  high-Q  photonic  cavity.  The  held  intensities  inside  the  structure  are  shown 
in  Fig.  A- 2b.  The  first  five  holes  next  to  the  cavity  are  quadratically  tapered  to 
ameliorate  the  mode  mismatch  between  the  cavity  mode  and  the  Bragg  mode  and  to 
increase  the  quality  factor  (Fig.  A-2c).  The  cavity  designs  are  optimized  numerically 
using  the  open-source  FDTD  package  MEEP  [127],  with  a  grid  spacing  of  a/20  and  a 
time  step  of  a/c.  We  normally  run  over  1,000  time  steps  until  the  mode  fully  decays. 
The  resonance  frequencies  and  quality  factors  are  obtained  through  a  filter  diagonal- 
ization  method  HARMINV  [128]  (error  5u/l o  <  10-3).  For  a  lattice  constant  of  a, 
the  optimized  cavity  design  has  a  width  w  =  2.8 a,  height  d  =  1.4a,  air  hole  length 
hx  =  0.52a  and  air  hole  width  hy  =  0.84ie.  We  obtain  a  cavity  with  a  resonance  at 
u  =  a/ X  =  0.447,  a  quality  factor  ( Q )  of  more  than  107,000  and  an  ultra-small  mode 
volume  of  1.36(A/n)3,  where  n  is  the  refractive  index  of  the  polymer. 

A. 2. 3  Theory 

To  determine  the  minimum  concentration  that  can  be  sensed  with  this  device,  we 
begin  by  considering  the  relationship  between  the  concentration  of  the  gas  and  the 
change  in  the  cavity  resonance:  ^  oc  AA  oc  C  [120],  where  A  A  is  the  wavelength 
shift,  A  is  the  cavity  resonance  wavelength,  AV  is  the  change  in  volume,  V  is  the 
total  volume  of  the  structure,  and  C  is  the  concentration  of  the  analyte.  Filling 
in  proportionality  constants,  this  can  be  written  as  ^  =  C\  A/-  =  C1C0C,  where 
Ci  quantifies  the  relationship  between  volume  change  and  wavelength  shift  and  Cq 
gives  the  relationship  between  concentration  and  swelling.  The  minimum  detectable 
resonance  shift  gives  the  minimum  detectable  concentration.  A  shift  in  the  resonance 
wavelength  can  be  measured  either  by  internal  luminescence  of  the  polymer  (using 
embedded  emitters)  or  by  a  passive  transmission  or  reflection  measurement.  Assuming 
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Figure  A-2:  (a)  Band  diagram  showing  four  distinct  modes  of  a  one-dimensional 
infinite  periodic  waveguide  comprised  of  a  unit  cell  with  a  period  of  a.  The  width  of 
the  polymer  beam  is  w  and  the  air  hole  has  dimension  hx  x  hy.  (b)  The  electric  held 
intensity  (|£'|2)  of  the  fundamental  mode  in  the  ID  nanobeam  cavity  in  the  z  =  0 
and  y  =  0  planes.  A  cross  section  of  the  Gaussian  intensity  profile  is  also  shown,  (c) 
Detail  of  the  lattice  constant  taper  used  to  create  a  high  Q  cavity. 


shot-noise  limited  detection, 

(j  .  _  i 
mm  “  qVtt 

where  £  =  1/coCi,  Q  is  the  quality  factor,  T  is  the  detected  photon  flux  (counts/s), 
and  r  is  the  integration  time  of  the  measurement. 

From  similar  work  with  polymer  coatings,  we  estimate  Co  ~  5  x  102  [120].  For 
volumetric  expansion,  c\  —  however  we  expect  that  at  low  concentrations,  the 
surface  of  the  structure  will  expand  more  than  the  interior,  causing  some  amount  of 
non-uniform  swelling  and  local  strain.  We  performed  FDTD  simulations  and  obtained 
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Ci  =  0.11  for  our  device.  Simulations  of  our  device  design  result  in  Q  values  above 
105.  For  this  calculation  we  use  T  «  1000  counts/s  on  a  single  photon  detector,  which 
was  typical  in  our  measurements.  This  yields  a  theoretical  concentration  sensitivity: 


StH  = 


qV f 


=  6  ppb/\/Hz. 


(A.2) 


A. 2. 4  Device  fabrication 

To  construct  the  devices,  we  spin-coat  a  10  nm  thick  layer  of  water-soluble  polyvinyl 
alcohol  onto  a  silicon  substrate,  followed  by  400  nm  of  PMMA.  The  patterns  were 
defined  using  EBL,  then  immersed  in  water,  allowing  the  PMMA  him  to  float  to 
the  surface,  where  it  could  be  transferred  to  any  rough  or  patterned  carrier.  Figure 
A-3  shows  two  examples  of  structures  onto  which  these  polymer  him  devices  can  be 
transferred,  including  tissue  consisting  of  minimally  absorbent  organic  hbers  with 
large  air  gaps  (Figs.  A-3a  and  A-3b).  These  were  used  to  lift  the  hlms  out  of  the 
water.  Arrays  of  devices  were  spaced  such  that  ~  80%  were  suspended  over  air  gaps 
when  the  him  is  randomly  oriented  on  the  tissue.  It  is  also  possible  to  suspend  the 
him  across  a  gap  formed  by  two  pieces  of  silicon  secured  to  a  glass  slide  (Fig.  A-3c). 
This  design  was  ideal  for  gas  sensing  experiments,  as  a  lens-corrected  objective  could 
be  used  to  image  through  the  glass  slide.  Figure  A-3d  shows  an  SEM  micrograph  of 
a  fabricated  device,  showing  the  smooth  edges  of  the  structure.  Additionally,  wire 
hoops,  plastic  transparency  paper  with  cut  holes,  rough  dense  paper  materials,  and 
tapered  hber  optic  lines  were  successfully  used  to  suspend  these  structures. 


A. 2. 5  Experimental  seutp 

We  tested  the  devices  on  a  confocal  microscope  setup  with  independent  pump  and 
collection  paths.  To  facilitate  the  optical  characterization  of  the  cavity  resonances, 
we  doped  the  PMMA  with  Coumarin  6  dye,  5%  by  weight.  The  dye  molecules  were 
excited  with  a  405  nm  laser  and  emit  from  450  nm  to  650  nm.  Figure  A-4a  presents  a 
cavity  Q  factor  in  excess  of  104.  Initial  tests  analyzed  IPA  vapor,  though  any  gas  or 
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Figure  A-3:  (a)  PMMA  film  on  spun  tissue  paper  (10 x  mag.)  (b)  Detail  of  structure 
in  PMMA  film  on  spun  tissue  paper  (lOOx  mag.)  (c)  Test  structure  using  two  pieces 
of  silicon  to  form  air  gap.  (d)  SEM  micrograph  of  fabricated  structure. 


vapor  that  interacts  with  the  polymer  film  can  be  used.  The  IPA  vapor  concentration 
in  the  system  was  controlled  by  a  dilution  gas  flow  setup.  The  sample  was  imaged 
through  a  glass  cover  slip,  using  a  lens-corrected  objective.  The  volume  containing 
the  gas  and  sample  was  roughly  10  mL.  A  gas  delivery  line  and  a  vent  line  allowed 
precise  control  over  the  gas  pressure.  Once  mounted  in  the  gas  cell,  the  polymer 
cavity  was  exposed  to  15  seem  dilute  gas  flow  for  5  minutes,  then  exposed  to  air  flow 
for  20  minutes  for  recovery. 


A.  2. 6  Discussion 

Four  gas  concentrations  were  tested,  and  the  resulting  shifts  are  detailed  in  Fig.  A-4b. 
The  three  lowest  concentrations  seem  to  follow  a  linear  trend,  while  the  highest  con¬ 
centration  appears  to  have  saturated  the  device.  When  a  gas  infiltrates  the  polymer 
cavity,  two  effects  occur:  the  refractive  index  changes  and  the  material  swells.  PMMA 
can  absorb  up  to  2wt%  of  its  mass  [129,  130].  The  refractive  index  of  the  tested  gas 
(niPA  =  1.38)  at  this  full  saturation  concentration  would  cause  a  ~  1  nm  blue-shift  in 
a  cavity  with  a  resonance  around  630  nm,  which  is  much  smaller  than  the  ~  22  nm 
redshift  caused  by  the  swelling  of  the  material  under  the  same  conditions,  making  the 
polymer  deformation  the  dominant  effect  [120].  We  show  full  recovery  of  the  sensor 
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Figure  A-4:  (a)  Photoluminescence  spectrum  of  a  higli-Q  cavity.  The  resonance  is 
at  627  nm  and  the  quality  factor  is  around  13,  000.  (b)  Resonance  shifts  plotted 
against  gas  concentration.  The  point  with  the  highest  concentration  seems  to  show 
saturation  behavior,  (c)  After  the  stimulant  gas  is  removed  from  the  system,  the 
resonance  returns  to  it’s  original  position.  The  baseline  signal  decreases  as  time 
progresses  as  a  result  of  bleaching  of  the  dye  in  the  PMMA. 


after  the  1PA  vapor  has  been  removed  from  the  system  (Fig.  A-4c).  Additionally, 
we  anticipate  that  low  concentrations  of  vapors  that  corrode  the  polymer  can  also  be 
detected,  such  as  ozone  and  acetone  vapor. 

When  the  cavity  is  exposed  to  the  minimum  concentration  tested  (600  ppm),  we 
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measure  a  0.6  nm  wavelength  redshift  from  631.9  nm  to  632.5  nm.  The  quality  factor 
of  this  cavity  was  2  x  103.  The  Q  may  be  improved  by  using  a  higher  resolution 
e-beam  lithography  tool  with  image-stitching  capabilities.  These  values  correspond 
to  an  experimental  sensitivity  of  Sexp  =  10  ppm/v^Hz-  From  Eq.  A. 2,  we  expect  a 
cavity  with  this  Q  factor  to  have  a  sensitivity  of  Sth  =  287  ppb/\/Hz,  which  is  over 
an  order  of  magnitude  better  than  the  performance  we  observed.  There  are  several 
possible  explanations  for  this  discrepancy:  first,  our  value  for  Ci  is  based  on  simulation 
and  assumes  an  even  expansion  of  the  surface  of  the  structure.  In  reality,  while 
every  effort  was  made  to  ensure  a  smooth,  even  gas  flow,  some  parts  of  the  structure 
may  have  been  oriented  such  that  the  top  surface  may  have  swelled  more  than  the 
bottom.  Second,  the  swelling  coefficient  Co  cannot  be  experimentally  measured  using 
our  setup  and  our  polymer  cavity  might  have  a  slightly  different  swelling  coefficient 
than  in  [120].  Given  the  precision  of  the  gas  delivery  and  vent  hoses,  extremely  low 
gas  concentrations  are  not  possible  to  verify  using  our  current  testbed.  Finally,  we 
expect  that  sources  of  noise  exist  in  our  experimental  setup  beyond  shot-noise,  such 
as  vibrational  and  thermal  effects.  We  also  see  a  decrease  in  our  cavity  Q  when  the 
gas  is  introduced  into  the  system.  This  can  be  accounted  for  by  the  increased  swelling 
of  the  surface  of  the  material;  surface  swelling  changes  the  relationship  between  a  and 
hx  and  between  w  and  hy,  as  detailed  in  the  inset  of  Figure  A-2a. 

In  conclusion,  compact,  low-cost,  liigh-Q  ID  polymer  photonic  crystals  with  high 
sensitivity  provide  the  opportunity  for  state-of-the-art  gas  sensing  on  almost  any 
rough  or  patterned  surface  imaginable.  We  have  demonstrated  the  implementation 
of  these  sensors  and  their  recovery  after  stimulant  gases  are  removed  from  the  sys¬ 
tem.  We  expect  that  experimental  sensitivity  will  increase  with  improved  fabrication 
processes  and  with  improved  readout  techniques,  such  as  doping  with  an  agent  that 
will  provide  stronger  fluorescence  or  with  the  implementation  of  cavity  transmission 
measurements.  This  platform  could  be  useful  in  a  variety  of  industrial,  environmen¬ 
tal,  and  even  biological  settings.  These  devices  have  the  potential  to  be  applied  to 
a  variety  of  sensing  applications,  especially  with  surface  functionalization  for  specific 
stimulants. 
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A. 3  High-Q  suspended  polymer  nanocavities  for  pre¬ 

cision  force  and  displacement  sensing 
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